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Abstract 

We consider the cubic defocusing nonlinear Schrodinger equation in the two dimensional 
torus. Fix s > 1. CoUiander, Keel, Staffilani, Tao and Takaoka proved in [CKS"*"!!] the 
existence of solutions with s-Sobolev norm growing in time. 

We establish the existence of solutions with polynomial time estimates. More exactly, 
there is c > such that for any ^ 1 we find a solution u and a time T such that 
— ^II'^(0)I|h-'- Moreover, time T satisfies polynomial bound < T < /C^. 
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1 Introduction 

Let us consider the periodic cubic defocusing nonlinear Schrodinger equation (NLS), 



{— idtu + An = |tipu 
u{0,x) = uoix) 

where x E = M7(27rZ)2, t G M and u : M x ^ C. 

The solutions of equation ([T]) conserve two quantities: the Hamiltonian 



(1) 



E[n]{t) 

and mass 



M[u]{t)= \u\'^dx{t)= \u\'^dx{0), (2) 

which is just the square of the L^-norm of the solution for any t > 0. It is useful to study 
solutions u{t) in a family of Sobolev spaces H'^ with the corresponding ff*-norms 



where (n) = (1 + |np)"^/^ and, 

u{t,n) := / n(t,x)e~*"'^ dx. 




The local-in-time well-posedness for any uq £ i?*(T^), s > was proven by Bourgain jBou93j . 
This along with the two conservation laws, implies existence of a smooth solution M for all time. 
It follows from conservation of energy E[u\{t) that the iif^-norm of any solution of ffl is uniformly 
bounded. Our main goal is to look for solutions whose higher Sobolev norms \\u{t)\\jjs(^j2^, s > 1, 
can grow in time. 

If the ff^'-norm can grow indefinitely for some given s > 1, while the ff^-norm stays bounded, 
then we have solutions which initially oscillate only on the scales comparable to the spatial period 
and eventually oscillate on arbitrarily small scales. To see that compare these norms. The only 
possibility for H'^ to grow indefinitely is that the energy of a solution of ([T]) can penetrate to 
higher and higher Fourier modes. 
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On the one-dimensional torus, equation ([T]) is completely integrable due to the famous result 
of Zakharov-Shabat |ZS71j (see also |GKP12) ). AsacoroUary ||^i(t)||Hs(Ti) < C'II^(0)IIh''(ti)) ^ ^ 
1 for all t > 0. If one replaces the nonlinearity \u\'^u = duP{\u\'^) in (l]J) with a more general 
polynomial, then Bourgain |Bou96| and Staffilani |Sta97a| proved at most polynomial growth of 
Sobolev norms. Namely, for some C > we have 

\\uit)\\H'' <t'^^'~'^^\\u{0)\\H'' for t^oo. 

In |BouOOaj Bourgain applied a version of Nekhoroshev theory. He proved that for a 1-dimensional 
NLS with a polynomial nonlinearity satisfying P{0) = P'{0) = P"{0) = for s large 

and a typical initial data u(0) G H'^{T) of small size e, i.e. ||m(0)|| < e we have 

sup < Ce, 

\t\<T 

where T < with A = A{s) — )• as s — )• +oo. This is an indication of absence of a polynomial 
growth and motivated Bourgain [BouOObj to pose the following question: 

Are there solutions in dimension 2 or higher with unbounded growth of H'^ -norm for s > 1? 

Moreover, he conjectured, that in case this is true, the growth should be subpolynomial in 
time, that is, 

< i^||ii(0)||Hs for t oo, for all e > 0. 

There are several papers obtaining improved polynomial upper bounds for the growth of Sobolev 
norms for equation ([T]) and also generalizing these results to other nonlinear Schrodinger equa- 
tions either on M, or M^, or on compact manifolds |Sta97bl ICiDKSnH IBoun41 IZhoOSI K^WlOl 
ISohlH ICK012] . Similar results have been obtained for the wave equation |Bou96j and for the 
Hartree equation [SohlObl ISohlOa] . 

All of the cited above papers give upper bounds of the growth but do not obtain orbits which 
undergo growth. Indeed, there are few results obtaining such orbits. In [Bou96j, Bourgain 
constructs orbits with unbounded growth of the Sobolev norms for the wave equation with a 
cubic nonlinearity but with a spectrally defined Laplacian. In |GG10llPocll] . it is shown growth 
of Sobolev norms for the Szego equation, and in [Pocl2] for certain nonlinear wave equation. 

Concerning the nonlinear Schrodinger equation, Kuksin in [Kuk97b] (see related works 
|Kuk95l IKuk96[ IKuk97a[ IKuk99j ) studied the growth of Sobolev norms but for the equation 

—iw = —6Aw + \w\'^^w, 6-^1, p > 1. 

He obtained solutions whose Sobolev norms grow by an inverse power of 6. Note that us{t, x) = 
5~2w{d~^t, x) is a solution of ([l]). Therefore, the solutions that he obtains correspond to orbits 
of equation ([T]) with large initial data. The present paper is closely related to |CKS"'"10] . In this 
paper, it was shown that for any s > 1 the H^-noicm can grow by any predetermined factor. 
The initial data there are not required to be large as |Kuk97b] . but rather have a small initial 
H^-noim with s > 1. Essentially using construction from this paper lOj we not only 

construct solutions with similar properties, but also estimate their speed of diffusion. 
The main result of this paper is 

Theorem 1. Let s > 1. Then there exists c > with the following property: for any large 
/C ^ 1 there exists a a global solution u{t, x) of ([Ip and a time T satisfying 

< T < /C^ 
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such that 

\\u{T)\\h^>1C\\u{0)\\hs. 
Moreover, this solution can be chosen to satisfy 

< ^-(s-l)c/4+2/(s-l)_ 



Note that Theorem [T] does not contradict Bourgain conjecture about the subpolynomial 
growth. Indeed, Theorem [T] only obtains solutions with arbitrarily large but finite growth in the 
Sobolev norms whereas Bourgain conjecture refers to unbounded growth. 

Remark 1.1. Even if Theorem^ is stated for ([T]) in the two torus, it can be applied to the d 
dimensional torus with d >2, since the solution we obtain is also a solution for equation ([T]) in 
the T'^ setting all the other harmonics to zero. 

Remark 1.2. In fact, we can obtain more detailed information about the distribution of the 
Sobolev norm of the solution u{T) from Theorem^ among its Fourier modes. More precisely, 
we can ensure that there exist ni,n2 G such that 

\HT)\\hs > |ni|2^|n„,(r)|2 + |n2p^|n„,(r)|2 > JC. 

That is, when t = T the Sobolev norm is essentially localized on two Fourier coefficients. 

Remark 1.3. Using more careful analysis of the proof we can establish existence of solutions 
whose Sobolev norms are lower bounded for each time t S [1,7"]. Namely, 

in\\uit)\\Hs > ^^+ln||n(0)||H- 

Remark 1.4. Our solutions differ from solutions studied in ]CKS^ id^ in a substantial way. If 
one applies to information about dynamics contained in \CKS^ IC^ supplied with the theory of 



normal forms and a beautiful trick of Shilnikov SiWV, then it is possible to compute certain 
"local maps" and diffusion time. It turns out to be super- exponential in fC, namely, it grows 
as C'^" for some C > and a > 2 (see Section 2.2 for more details). Even equipped with 



the aforementioned dynamical technique in order to obtain polynomial diffusion time we need to 



achieve ~ In /C cancelations. These cancellations are spilled out in Section 2^ on an heuristic 
level and then worked out in Sections\^ and\^ 

In |CKS"'"10] initial conditions of solutions with growth of Sobolev norms can be chosen with 
small ||M(0)||jys|^ In our case it is also possible, but leads to slowing down of time of growth. 
This fact is explained in Appendix [C| 

The present paper deals with growth of Sobolev norms for a Hamiltonian partial differential 
equation. We show the existence of unstable solutions. As we have explained, there have not 
been many results showing the existence of these instabilities. In |CE10] a solution of ([T]) 
with spreading of mass among modes is constructed. Nevertheless the spreading does not lead 
to growth of Sobolev norms. In [Hanllj a progress toward infinite growth of Sobolev norms is 
made. Let us say also that in the past decades there has been a considerable progress in the study 
of other types of dynamics for Hamiltonian partial differential equations. For instance, in the 
existence of periodic, quasi-periodic or almost-periodic solutions (see e.g. |Rab78t |Way90[ W931 



'As Terence Tao pointed to us, our solutions have small L^-norm, but not //'-norm 
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IKPn.Sl IKuk93[ IKP961 IBerOTl IBBllj ^. in Nekhoroshev type results (see e.g. |Bam97[ IBam99p 
and normal forms (see e.g. |Bam03l IBG06L IG1P09L IGKP121 IPP12j '). Of particular interest for 
the present paper are |Bou98t lEKlO] since, in these papers, the authors study the existence 
of quasi-periodic solutions for the nonlinear Schrodinger equation in the 2-dimensional torus 
|Bou98j and in a torus of any dimension |EK10j . Nevertheless, they consider slightly different 
equations containing a convolution potential. 



2 Main ideas and structure of the proof 

One of remarkable contributions in [CKS^IO] is the formulation of a finite-dimensional toy 
model, which after a certain lift approximates solutions of ([T]). The Hamiltonian of the toy model 
from jCKS"'"10] has a specific form. It has a nearest neighbors interaction and is integrable inside 
a certain family of 4-dimensional planes. In this section we present a class of Hamiltonians with 
a nearest neighbors interaction for which our method applies. It is specified at the end of Section 



2.1 Features of the model 

• Write ([T]) as infinitely ODE's for Fourier coefficients of solutions. It is a Hamiltonian 
system with Hamiltonian T-L (see Q). 



{Two step reduction) 

— Obtain a Normal Form of the original Hamiltonian near the origin by removing non- 
resonant terms (see Theorem [2]). 

— Use gauge freedom to remove linear and some non-linear terms (see (|13|)). 



{The Toy Model) 

Select a finite subset of Fourier coefficients A in 1? so that they can be split into pairwise 
disjoint generations A = U^^^Aj and only neighboring generations Kj and Aj+i interact. 

This can be done so that dynamics of each element in each generation has exactly the 



same as dynamics of any other member of this generation (see Corollary 3.2). Truncating 



we are reduced to a complex A^-dimensional system given by a Hamiltonian 

N N-1 

where each hj is complex valued, and the symplectic form = ^dbj A bj. The system 
conserves mass A4{b) = J2f=i l^iP- We study the dynamics restricted to mass A4{b) = 1. 
Dynamics of this Hamiltonian is called in [CKS^IO] the Toy Model and is the focal point 
of analysis. It is convenient to study this system in real coordinates and identify C = M^. 

Notice also that the Hamiltonian h{b) can be viewed as a Hamiltonian on a lattice Z with 
nearest neighbor interactions. Our main result relies on the construction of energy transfer 
from 63 ~ 1, bj ~ 0, j 7^ 3 to b]\f^2 ^ ^, i 7^ ~ 1 for this Hamiltonian. Construction of 
a somewhat similar energy transfer for the pendulum lattice is done in [KLSllj . 
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Figure 1: Planes approximating solutions 

• {Invariant low-dimensional subspaces) 
Notice that each 4-dimensional plane 

Lj = {bi = ■ ■ ■ = = bj+2 = ■ ■ ■ = bN = 0} 

is invariant. Moreover, dynamics in Lj is given by a simple Hamiltonian 

Denote Aij{bj,bj+i) = \bj\'^ + Both hj and Aij are conserved. The mass Aij is 

assumed to be 1. 

The solutions constructed stays close to the planes 
and go from one intersection Ij = Lj n ij+i to the next one 
Ij^i = Lj^i n Lj^2 consequently for j = 3, . . . , N — 2 (see Figure^. 

To make a closer look at solutions we need to understand dynamics in the planes Lj^s. 

• {Integrable dynamics in each plane Lj) 

Dynamics in each 2-dimensional plane Lj is integrable. Indeed, there are two first integrals 
hj and Mj in involution. By Arnold-Liouville theorem away from degeneracies the 4- 
dimensional plane Lj is foliated by 2-dimensional invariant tori with dynamics smoothly 
conjugated to a constant flow. 

We are interested in two specific periodic orbits: 0j-direction {\bj\ = 1, = 0} and 

0j+i-direction = 1, 6j = 0} and in a family of heteroclinic orbits {7j} connecting 

the former with the later. All these orbits can be found explicitly, but their existence can 
be predicted having hj and Mj satisfying some properties. 

— Having the mass Mj = \bj\'^ + l^j+iP conserved it is natural to expect that the 
boundary is invariant. The boundary consists of bj = and 6j+i = (both periodic 
orbits) and belong to the same /ij-energy surface. 
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— It is a straightforward calculation to check that both orbits are hyperbolic, i.e. of 
saddle type. 

— Notice that {hj = |, A4j = 1} is a 2-dimensional surface with the boundary given 
by periodic orbits bj = and = 0. Away from these periodic orbits it is a locally 
analytic surface, i.e. gradients Vhj and V-Mj are linearly independent. 

— Away from the periodic orbits bj = and bj+i = the surface {hj = \, Mj = 1} 
consists of stable and unstable 2-dimensional manifolds. Unless the periodic orbits 

bj = and = on {hj = |, ^Aj = 1} arc separated by a degenerate periodic 
orbit, they have to be connected by these manifolds. 

— Now we verify that there is no such a degenerate periodic orbit. Moreover, we find 
explicitly the family of connecting heteroclinic orbits. Even though these explicit 
formulas is not used in our proof. 

Write in polar coordinates bk = \/T^e^^^^ k = j,j + 1. The mass conservation becomes 
^Aj{b) = Tj + Tj+i, the symplectic form Q = -^drj A dOj and the Hamiltonian 

hj ( e^'^ , e^'^+i) = ^ [r| + r|+i + 4r,T,-+i cos 2(% - . 

Then the equation of motion are 

Oj = Vj — 2rj+i cos 2{9j — 9j+i) 
9j-^-l = rj+i — 2rj cos2{6j — ^j+i) 

fj = Avj rj^i sin 2{6j — 6j+i) 
■Tj+i = -4r-j Tj+i sm2{9j - 9j+i). 

For the energy surface hj = ^ we have 

— Two families of periodic solutions 

{{ej,ej+i,rj,rj+i) : Vj = 0} and {(6*^, ^j+i, r^, r^+i) : r^+i = 0}. 

— Each family has two special solutions: 2{6j — Oj+i) equals either ^ and Both 
planes are invariant: ^{Oj — Oj+i) = —{vj + rj+i)(l + 2cos2(0j — = 0. Denote 
Tj = {2(% - Oj+i) = ^ ( mod 27r), Vj = 0}. 

— On M.j = 1, hj = \, 9j — 9j+i = ^ we have Vj = VjVjj^i = —fj+i. Thus, there is a 
heteroclinic orbit 7j connecting Tj with the second family r^+i = 0. 

Now we can be more specific in location of orbits: 

The solutions constructed go from one periodic orbit T2 to the next T3 
along 72, then from T3 to T4 along 73 and so on for j = A, . . . , N — 2. 

In a view of the above discussion we have the following description: 

e,^o,\i-j\>i ^,«o, > 1 (4) 
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Figure 2: Heteroclinic orbits 



{Local behavior of periodic orbits Tj) Due to the above analysis, the periodic orbits Tj 
viewed in have at least two expanding and two contracting directions: one pair 
from Lj_i-plane and the other from Lj-pane. Due to symmetry of the restricted systems 
in Lj_i-plane and Lj-plane these periodic orbits have multiple hyperbolic eigenvalues. 
Multiplicity turns out to be exactly 2. 

{Resonant normal forms nearTj) Presence of resonance complicated analysis and as for- 



mulas (67) show resonance changes loc al beh avior compare to the linear case. To resolve 
it we use a beautiful trick of Shilnikov and obtain precise information about local 



behavior, which is explained in Section 2.2 



• ( Connecting heteroclinic orbits) As we showed above there are orbits 7^ connecting Tj with 
Tj+i for each j = 3, . . . , n — 2. We need to analyze dynamics near these heteroclinic orbits. 

• (Local almost product structure) Once we obtain information about behavior near Tj's and 
near connecting orbits jj , we can describe dynamics of the Toy Model as if it close to the 
direct product of (N — 3) planes Lj, j = 3, . . . , N — 1. 

Properties of the Hamiltonian h{b) used in the proof. 

As we mentioned in the introduction to this section we do not use a specific form of h. Here is 
the list of properties that we need. 

• h has nearest neighbors interaction; 

• h has 2-dimensional (complex) invariant planes intersecting transversally; 

• there are two first integrals (coming from two conserved quantities: energy and mass); 

• some generic properties of h and Ai. 
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2.2 The dynamics close to the periodic orbits: a heuristic model 

One of the crucial steps in analyzing the toy model h{b) is the study of the dynamics in a 
neighborhood of the periodic orbits Tj. Namely, we want to analyze how points which lie 
close their stable invariant manifold evolve under the flow until reaching points close to their 
unstable one (see Figure |3]). As we have explained, these periodic orbits are of mixed type 
(four eigenvalues are hyperbolic and the rest are elliptic). Since in each plane Lj dynamics is 
the same explained in the previous section, the hyperbolic eigenvalues have multiplicity two 
and, therefore, are equal to A, A, —A, —A for some A > 0. Since in this section serves exposition 
purposes we let A = 1 and set the elliptic modes to zero. |^ 
Essentially the study has three steps: 

• Using conservation of Ai, make a simplectic reduction so the periodic orbit becomes a 
fixed point. 

• Perform a normal form procedure to reduce the size of the higher order non-resonant terms. 

• Analyze the dynamics of the new vector field and achieve a cancelation for a local map. 



The first step is performed in Section 4.1 It leads to a Hamiltonian of two degrees of freedom 
of the form 

H{p,q) =Piqi+P2q2 + H4{q,p), 

where H4 is a homogeneous polynomial of degree four. The variables {pi,qi) correspond to the 
variable bj-i after diagonalizing the saddle and the variables (^2,^2) correspond to bj+i. 

Fix a small a > 0. To study the local dynamics, it suffices to analyze a map from a section 
S- = {qi = a, \pi\, \q2\, \p2\ ^ o"}, to a section 11+ = {p2 = a, \pi\, \qi\, \q2\ ^ cr} (see Figure 
[3]). Using rescaling assume a = 1. This can change time by a fixed factor. 

Since we are in a neighborhood of the origin, one would expect that the dynamics of the 
system associated to this Hamiltonian is well approximated by its first order, that is, by a linear 
equation. Then, the solutions are just given by 

pi(t) =p?e*, ^i(t) = gOe-* 
P2{t)=ple\ q2{t)=qle-' 

and then the local map Bq from U C 5]_ to for this system sends points 

(P?,9?,P2,'?2)~(5,1,AV5) 

to 



where < 5 ^ 1. Moreover, the travel time of orbits by this map is always T = - In 0(1). 



*To be more precise near each saddle, the elhptic directions remain almost constant and, since they will be 
taken small enough, it turns out they do not make much influence in the dynamics of hyperbolic components. 
Thus, to simplify the exposition, we set the elliptic modes to zero and study how the hyperbolic ones evolve. This 
implies that we only need to study three modes hj and 6j+i. This analysis is performed in Section[5]in great 

detail. 
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Figure 3: Local map 



We will see that the image point changes substantially when we add to the system, due to 
both resonant and nonresonant terms. To exemplify this, we consider a simplified model which 
in fact contains all the difficulties that the true model has, 

H{p, q) =Viqi+ V2q2 + qlvl + PipI- (5) 

Since the term p\p2 is nonresonant, we first perform one step of normal form (x, y) = '^{p,q) 
(see Section [5] for details). It can be easily seen that the change ^ is of the form 

q) = {pi, qi + 0{pipl),p2, q2 + 0(p?P2)) (6) 
and, therefore, keeps the size of initial points of the form 

(p?,g?,P2,^2)~(<5,l,V5,V5). 

That is, = '^{p^,q°) satisfies 

{x^^ylx^y'^) - {s,l,V6,V6y 

The change to normal form leads to a Hamiltonian system of the form 

H'{x, y) = xiqi + X2y2 + yix\ + higher order terms. 

Drop the higher order terms. Then, the solutions of the system associated to this Hamiltonian 
can be computed explicitly and are given by 



Xl = 




+ 2y',{xl)He' = 


(x? + 2y?(x°)2: 


t) e* 


yi = 




■t 






X2 = 










y2 = 


yy 


* - 2(y?)^x°ie*. 
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Thus, since the travel time is t = — In ^fS + 0(1), it is clear that the nonlinear terms are bigger 
than the linear ones, leading to an image point of the form 



f f f f 

X It If* 1 



\/51n(l/(5),\/5,l,51n(l/5)) . 



Using ([g]), in the original variables the image point of the map Bi associated to Hamiltonian H 
is of the form 

Bi (p?, qlvl q°2) ~ ln(l/5), Vd, 1, 6 ln\l/6)) . 



We want to emphasize that the presence of these logarithmic terms is a serious problem we need 
to deal with. Recall that we need to travel through — 3 saddles (T3 T4 . . . TAr_i). 
Roughly speaking, this implies that we need to compose — 4 local maps. Thanks to the 
symmetries, at each saddle we can consider a system of coordinates such that the dynamics is 
essentially given by a Hamiltonian of the form ([5]). Moreover, since at each local map we gain 
some logarithms, the initial points of the local map associated to the j saddle are of the form 

which, thanks to ([6]), in the normal form variables satisfy 

Then, proceeding as before, these points are mapped to points of the form 



f f f f 

X 01 1 

1 ) j/l 5 -''2 ' ^2 



which in the original variables read 

Bi [plqlp^ql) ~ (^5\n^'-\l/5),V~5A,5\n^\l/5)) . 



That is, the amount of logarithms doubles at each step and thus grows exponentially. This 
accumulation of logarithmic terms leads to very bad estimates. Indeed, to keep track of the 
orbit after — 4 local maps, we would need that 

51n2^"'(l/5) < 1. 

Therefore, we would need to choose 5 extremely small with respect to N . 

For example, if 5 > C~'^ " ~ C"^" for some C > independent of A^, then the above 
expression gives 

C''^^\r'' lnC7)2"'%>lfora<l. 

In this case, the constant A appearing in Theorem [4] would need to satisfy A ^ for some 6 > 
and independent of A^. As a result. Theorem 3 gives a diffusion time T ~ X'^K-fN In 1/5 > C'^ 
(see formula (22)). Thus, choosing such a small 6 would lead to very bad estimates for the 



diffusion time of Sobolev norms as we pointed out in Remark 1.4 



To overcome this problem, we modify slightly the initial conditions. Notice that if we choose 

x?-2y?(x^)^ln \/5 = 0, 



Xi such that 
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we obtain that at the end x{ ~ and thus we avoid the logarithmic term. This cancelation will 
be crucial in our proof. If we restrict Xj* to this set, we are taking (51n(l/(5) and therefore 
we will be sending points 

to points 

[x(,y(,xi,yl) ~ (V5,V6,l,Sln{l/6)) . 



The map will keep the same form expressed in the original variables, and, therefore, we will 
avoid having increasing separation from the invariant manifolds. 

2.3 Outline of the Proof 

• Find symplectic coordinates near the origin in where the original Hamiltonian T-l sim- 
plifies (see Theorem [2|. Namely, T-L oF = T? + Q + TI, where P is a quadratic Hamiltonian, 
G is of degree four and only contains resonant terms, and TZ is smaller. 

• Dynamics of D + G has invariant finite-dimensional subspaces, which give rise to a sim- 



pler (no simple!) finite-dimensional Hamiltonian h{b) given by (19). In terminology of 
|CKS"'"10] this Hamiltonian defines the Toy Model. In Theorem [3 we obtain orbits of the 
toy model which have transfer of energy. 

We show that are solutions of the system associated to Ti which are close to those of the 
toy model for long enough time (Theorem |4]) . These orbits undergo the wanted growth of 
the Sobolev norm. 

The proof of Theorem [3] occupies most of the paper. Theorems [2] and |4] are deferred to 
Appendices |A] and [B] respectively. Now we describe the plan of the proof of Theorem [3] 



Following |CKS+ lO] we detect a collection of periodic orbits of h(b), defined in 

(32), and heteroclinic orbits connecting them (see (33)). 



The whole proof consists in a careful analysis of dynamics near the union of these periodic 
orbits and their connecting orbits. Our analysis naturally splits into 

— local dynamics near periodic orbits {Tj}j^^ and 

— global dynamics near heteroclinic orbits {7^}^^^^. 

More formally, Theorem [3] follows from Theorem [5] The latter Theorem in turn follows 
from Lemmas 14.71 and 14. 



The Local Lemma 4.7 provides refined information about local behavior near periodic 



orbits {Tj}j with quantitative estimates. 

Global Lemma |4.8| provides refined information about local behavior near heteroclinic 



orbits from (33) with quantitative estimates. 



The proof of Local Lemma |4.7| consists of several steps. As we have explained in Section 
the periodic orbits {Tjlj have mixed type. Namely, in some directions the local 



2.1 



behavior is hyperbolic, while in others it is elliptic. It turns out that the closer orbits 
under investigation pass to the periodic orbits {Tj}j, the more decoupled (direct product- 
like) behavior they have. 
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In Section [5] we set all the elliptic variables zero and study the (4-dimensional) hyperbolic 
Toy Model. 

In Section [6] we use these results to deal with the full hyperbolic-elliptic system and prove 
Lemma 14.71 



• In Section [7] we prove Global Lemma 4.8 As we pointed out, this implies Theorem [5} 
which in turn, implies Theorem [Sj 

• Combining this result with Theorem [2| proved in Appendix [Aj and Theorem |4] proved in 
Appendix [B| we complete the proof of the main result (Theorem [T]) . 

We summarize this in the following diagram: 



Theorem [T] 

1t 



Theorem [2] 


+ 


Theorem [3] 


+ 


Theorem |4] 



1t 

Theorem [5] 

it 



(7) 



Local Lemma 14.71 



+ 



Global Lemma 14.81 



2.4 Major ingredients of the proof 

We summarize here the new set of tools that we apply to the problem compared to |CKS"'"10 

• In Theorem [2| we use a standard normal form (e.g. see |KP96j ). 

• Theorem [3] requires several new ideas: 

— Finitely smooth resonant normal form for hyperbolic saddles |BK94j . 



Shilnikov boundary value problem [Sil67| to study the local behavior close to the 
periodic orbits Tj. 



— As we explained for the model case in Section 2.2, to control the dynamics of the toy 
model we need a peculiar cancellation (see Section [s]) . 

— To have cancellations at each stage, we need to establish local product structure for 



the orbits we are interested in (see Definition 4.3). 



Due to the good control of the solutions of the toy model, we are able to approximate the 
solutions of the original systems with the ones of the toy model for longer time compared 
with jCKS+10] (see Theorem [i]). To achieve this, we also modify the set A (see condition 
6a). This modification allows to slow down spreading outside A. 
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3 The three key theorems 

We start the proof analyzing the infinite system of equations which describe the behavior of 
Fourier coefficients. Namely, consider the Fourier series of u, 

u{t,x) = a„(t)e^"'', Unit) := u{t,n). 

Therefore, the equation ([T]) becomes an infinite system of equations for {an}n&^j which are 
given by 



ni-n2+n3=n 

Note that this equation is Hamiltonian. Indeed, it can be written as 

dn = 2id n{a,a), 

where 

'H{a,a) =V{a,a) + g{a,a) (9) 

with 



1 2 1 1 2 
2 . ."-I lO-nl 



flni '^ras '^n4 • 



n\—n2+n3=nA 



We will study equation ([8| in a family of Banach spaces: all -ff^-Sobolev spaces with s > 1 as 
well as in the £^ -space. The £^ space is defined as 



a : — 7- C : ||a||^i = < oo 



Note that, £^ is a Banach algebra with respect to the convolution product. Namely, if a, 6 G £^ 
its convolution product a * 6, which is defined by 



{a*h)n= ^ ani^raz 
n\+n2=n 

satisfies 

||a* 5||^i < ||a||^i||6||^i. 

Finally, let us point out that the L^-norm conservation of ([T]), becomes now conservation of 
the £^-norm of a, defined as above. Namely, we have that ||a(t)||^2 = ||a(0)||£2 for all t G M. 

We want to study the evolution of certain solutions of equation ([s]), which will be small in 
the (.^ norm. Now we make an outline of the proof. 

The first step is to find out which terms make the biggest contribution to this evolution. To 
this end, we perform one step of normal form and bound the remainder in the £^-norm. 
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Theorem 2. For the Hamiltonian Ti in ^ there exists a symplectic change of coordinates 
a = r(a) in a neighborhood of in which takes it into its Birkhoff normal form up to order 
four, that is, 

-HoT = v + g + n, 

where Q only contains resonant terms, namely 

^(q!, tt) — — ^ ^ OLYi^C^ri2^ni^ni 

n\—n2+n2=nA 
|ni|2-|n2|2+|n3|2=|n4p 

and X-ji, the vector field associated to the Hamiltonian IZ, satisfies 

\\XnU<0{\\a\\%). 

Moreover, the change V satisfies 

||r-Id||^i <0{\\a\\%) . 

The proof of this theorem is postponed to Appendix |A| 

Once we perform one step of normal form, we have a new vector field 

- ian = \n\^an+ ^ an^a^anz + (hj^, (10) 

(ni,n2,n3)eylo(n) 

where 



I |2 I |2 I I |2 I |2l 

\ni\ - \n2\ + \n^\ = \n\ >. 



n, 



:iii 



As a first step, we focus our attention to the degree 4 truncation of it, which will give the main 
contribution to the dynamics. Namely, we consider the Hamiltonian 

u' = v + g, 

which has associated equations 

-id„ = |npa„+ ^ a„^a^a„3. (12) 

(ni,n.2,n3)e^o(n) 



Note that the £^-norm of a is a first integral of this system as well as for Q and (10). 
Namely, 

||a(t)||^2 = ||a(0)||^2 for ah t G M. 

Then, to study the dynamics of a close to the origin (in the £^-norm) we remove its linear 
terms using the variation of constants formula. Moreover, following [CKS"'"10) . we also remove 
certain cubic terms using the gauge freedom of equation ([T]) . To this end, we make the change 
of coordinates 

an = /3„e<^+l"l')*, (13) 
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where G > is a constant to be determined. The equations for (3 read 

^^/3n = GPn + ^ ^ /3ni /5n2 /^ns • 

{n\,n2,nz)'^Ao{n) 

Choosing G properly we can remove certain terms in the sum. Indeed, we spht the sum as 



- E 

(ni,n2,n3)eylo{»i) (»^l,»^2,'^3)e^o(»^) {ni,n2,n-i)&Ao{n) (ni,7i2,n3)Gv4o(n) {ni,n2,nz)&Ao(n) 
ni,n3^n ni=n n3=n ^1=713=71 



The last sum is just one term, which is given by — /3„|/3„p. The second and third sums, are in 
fact single sums and each of them is given by 



Recall that both (12) and (13) preserve the ^^-norm. Therefore, taking G = —2\\a\\'j2 = — 2||/3||^2 5 
we can remove these two terms. Thus, with this choice, we obtain the equation for (3, which 
reads 

-i^n = -/3„|/3n|'+ Yl /5niA^/3n3 (14) 

ni,n2,n3£A(n) 

where 

-^(n) = I ('^i; "•2,^3) G (Z^)^ : ni - 712 + n3 = n 

I |2 I |2 I I |2 I |2 / /I 

Pi| - P2I + PsI = \n\ ,ni n,ns n>. 
We define also the set of all resonant frequencies as 

A = |(ni,n2,n3,n4) G (Z^)^ : (ni, 712, 713) G ^(n4)|. 

Note that if 722, 71-3, 71,4) G A, then the four points form a rectangle in I?. 

We reduce this system to a finite-dimensional one, which corresponds to an invariant finite- 
dimensional plane. To this end, we consider a set A C Z^ such that the corresponding harmonics 
do not interact to the harmonics outside of A. Moreover, we obtain a set A such that the 
harmonics in A interact in a very particular way. This set was constructed in [CKSj lO]. We 
explain now its construction and impose an additional condition on A from [CKS^IO] . 

Fix > 1. Following |CKS+ 10] we define a set A C 7? consisting of N pairwise disjoint 
generations: 

A = AiU...UAAr. 

Define a nuclear family to be a rectangle {ni,n2,n3,n4) G A, such that rii and 713 (known as 
the parents) belong to a generation A-,- and 712 and 71,4 (known as the children) live in the next 
generation Aj+i. Note that if (^i, 7Z2, ?7-3, 71-4) is a nuclear family, then so are (7x1,77,4,77,3,77,2), 
(?^3, 7X2, 77,1, 77,4) and (7x3,7x4,771,772). These families are called trivial permutations of the family 
(7x1,7^2,^3,77,4). 

The conditions to impose to the set A are 

1a Closure If 7x1,7x2,7x3 G A and (7x1,7x2,773) G .4(7x), then rx G A. In other words, if three 
vertices of a rectangle are in A so is the last fourth one. 
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2a Existence and uniqueness of spouse and children For any 1 < j < and any ni G Aj, there 
exists a unique nuclear family (ni, 722, na, 714) (up to trivial permutations) such that ni is 
a parent of this family. In particular, each m G Aj has a unique spouse ns G Aj and has 
two unique children 712,724 G Aj+i (up to permutation). 

3a Existence and uniqueness of sibling and parents For any 1 < j < N and any n2 G Aj+i, 
there exists a unique nuclear family (?i.i, 77,2, 77,3, 714) (up to trivial permutations) such that 
ri2 is a child of this family. In particular each 712 G Aj+i has a unique sibling 714 G Aj_|_i 
and two unique parents 711,713 G Aj (up to permutation). 

4a Nondegeneracy The sibling of a frequency ti never equal to its spouse. 

5a Faithfulness Apart from the nuclear families, A does not contain any other rectangle. 

These are the conditions imposed on A in jCKS^lo] . We will impose an additional condition: 

6a No spreading condition Let us consider n ^ A. Then, n is vertex of at most two rectangles 
having two vertices in A and two vertices out of A. 

Proposition 3.1. Let /C ^ 1. Then, there exists N ^ I large and a set A C 1/^, with 

A = AiU.-.UAat, 

which satisfies conditions 1a - 6a and also 

Moreover, given any R > (which may depend on IC), we can ensure that each generation Aj 
has 2^"-^ disjoint frequencies n satisfying \n\ > R. 

The proof of Proposition 2.1 from [ CKS"'"10] applies. We prove a quantitative version of this 
proposition in Appendix [Cj 

We use the set A to obtain a finite dimensional dynamical system (of high dimension) ap- 



proximating (14). To this end, let us first note that, by Property 1a, the manifold 

M = 1^ G C^' : /3„ = for all 71 a| 



is invariant by the flow associated to (14) and is finite dimensional. Indeed, by Proposition 3.1 
its dimension is N2^~^. Equation (14) restricted to M reads as follows. For each 7Z G A we 



have 

|2 



if^n — /3n\(3n\ + '^f^n^hild-^ /^"childj /^"spouse + ^/^n parent -j^ /^^parcnt2 ng 

(16) 

Indeed, presence of parents, children, and the sibling are guaranteed by 2a and 3a- Note, that 
in the first and last generations, the parents and children are set to zero respectively. In fact, 
M has a submanifold of considerably lower dimension which is also invariant. 

Corollary 3.2. (cf fCKS^ 1(^ } Consider the subspace 

M = {/? G M : /3„j = f3n2 for all ni,n2 G Aj for some j} , 
where all the members of a generation take the same value. Then, M is invariant under the flow 



associated to (16) 
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The dimension of M is equal to the number of generations, namely N . To define equation 
(16) restricted to M, let us define 



bj = Pn for any n G Aj 



Then, (16) restricted to M becomes 

bj = -ih% + 2ibj {b]_, + b]^,) , j = 0,...N, 
which is a Hamiltonian system with respect to the Hamiltonian 



(17) 



(18) 



(19) 



J 3 

and the symplectic form Vt = ^dbj A dbj. 

Theorem 3. Fix a large 7^1. Then for any large enough N and 6 = e~"'^ , there exists an 
orbit of system (18), u > Q and Tq > such that 

and 



\bm\>^-5'' 



|6jv-i(To)| yl-S" 

\b,{n)\<5'' forj^N-1. 



Moreover, there exists a constant K > independent of N such that Tq satisfies 



0<To<KNln 



Remark 3.3. An analog of this proposition also holds for some smaller 5, e.g. 5 = C~ 
is related to Remark\1.4\ about time of diffusion without cancelations. 



(20) 
This 

Using ([17]), Theorem [I gives an orbit for equation (|14l). Moreover, both equations (|14l) and 

(21) 



(18) are invariant under certain rescaling. Indeed if b{t) is a solution of (18) 

b\t) = x-^{x-H) 



is a solution of the same equation. By Theorem |3] duration of this solution in time is 



T = X^Tq < X'^KjN'^, 



(22) 



where Tq is the time obtained in Theorem [3| which satisfies (20). 

We will see that, modulo a rotation of the modes (see (13)), there is a solution of equation 
( 10 ) which is close to the orbit of (|14l) defined as 



(23) 



I3^{t) = X'^bj [X~H) for each n G Aj 
I3^(t) = for each n A. 



To have the original system being well approximated by the truncated system, we need that 
A is large enough. Then the cubic terms in (10) dominate over the quintic ones. Nevertheless, 
the bigger A is, the slower the instability time by (22). Thus, we look for the smallest A (with 
respect to A^) for which the following approximation theorem applies. 
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Theorem 4. Let a{t) = {oin{t)}n&? solution of (10), li^{t) = 

solution of (14-) given by (23) and T be the time defined in (22). Suppose suppa(O) C A and 

a(0) = /3'^(0). Then, there exist a constant k > independent of N and 7 such that, for 



A = e 



K-yN 



and < t < T we have 



where G = -2||a(0)||22. 



\an{t) - e 



i{G+\n\^)t o\ 



(24) 
(25) 



Using the three key theorems: Theorems [2| [3] and |4] we complete the proof of Theorem [T} 

Proof of Theorem^ Using the change of variables F obtained in Theorem [2| from the solution 
Q obtained in Theorem |4] we define a = T{a), which is a solution of system (|8]). We show that 
this orbit has the properties stated in Theorem [T] 

To compute the growth of Sobolev norm of this orbit a, we use the notation 



Sj = |n|2^ for J = l,...,Af-l. 



(26) 



neA,- 



To estimate the mass of our solution recall that 2^ = 1 = |Aj|. We want to prove that 

\HT)\\h^^^ 



and estimate the mass ||a(0)||j;^2 of the solution. To this end, we start by bounding ||a(T)||/fs in 
terms of Sn~i- Since 



|a(T)||^. > V |n|2^|a„(T)p>57v-i inf |a„(r)p 

neA]v_i 



it is enough to obtain a lower bound for |ari(r)| with n £ A^v-i- Using the results of Theorems 
[2]and|4| we obtain 



|a„(T)| > |a„(r)| - |r„(a)(r) - an{T)\ 



> 



i{\n\^+G)T 



a„(T)-/3„^(r)e<W^+G)T 



(27) 



|r„(a)(r)-a„(r)|. 



We need to obtain a lower bound for the first term of the right hand side and upper bounds for 
the second and third ones. Indeed, using the definition of in (23) and the results in Theorem 
[3] we have that for n G Ajv_i, 

/3^(T)' = A-2|6jv-i(ro)p>^A-2, 



(the relation between T and Tq is established in (22)). 

For the second term in the right hand side of (|27|), it is enough to use Theorem |4] to obtain. 



a„(r)-/3„^(r)e<w^+^)^f < |a„(r)-/3„^(r)e<w^+^)^|) 



k-2 



< 
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For the lower bound of the third term, we use the bound for F — Id given in Theorem [2j Then, 

A-2 



|F„(a)(r) - a„(r)|2 < ||F(a) - a\\l < 



8 



Thus, we can conclude that 



A_-2 
2 



Now we prove that 



\am\ls<X-'S3 and ||a(0)||i, < 2^. 



(28) 



(29) 



By the definition of A in ( |24[ ), the second inequality implies that the mass of a(0) is small. On 
the contrary, the first inequality does not imply that the ff*-norm of a(0) is small. As a matter 
of fact ^3 is larg^ 

To prove the first inequality of ( |29[ ), let us point out that 



H- 



< H^" + (rn(a(0) - an(0))p 



We first bound ||q(0)||^s. To this end, let us recall that supp a = A. Then, 



a 



Using Theorem [4j we have that 



n 



2s 



/5n(0) + /3^(0) - an{0) 



Recalling the definition of (3^ in ( 23 ) and the results in Theorem [sj 

ElH''|/^n(0)f < (1-0^3 + ^'^E^^- 

neA 



S3 



From Proposition 3.1 we know that j 7^ 3, 



(30) 



Therefore, to bound these terms we use the definition of 5 from Theorem [s] taking 7 = 7(5 — 1). 
Since s — l>so — l>Ois fixed, we can choose such 7 » 1. Then, we have that 



neA 



< A-2S3. 



^As pointed out to us by Terence Tao. 
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From this statement, (30) and Theorem |4| we can conclude that 

\Ho)fHs<x-'S3. 

To complete the proof of statement ( |29| ) recall that the support of T(a) — a is 

= |n G : n = ni — 71,2 + n^, ni,n2,Ti^ G A| 



and apply Theorem [2j 

Using inequalities ( |28[ ) and (29), we have that 



HT)\\ 



HO)\\ls ~ ^3 



and then, applying Proposition |3.1[ we obtain 



ll«(0)|II,. ~ 2- 
It is left to estimate diffusion time T. Use Proposition 



3.1 



to set /C ~ 2(^^-1)^/2 j^j^jj ^ 



4:Hj/{s — 1), definition (24) to set A = e'^'^^ ~ /C'^/^^^"^). For time of diffusion we obtain 



|r| < K 7 X^N^ < IC 7 x:^/in2 41n2/C ^ 
~ In^ 2 (s - 1)2 - 



for large IC. This completes the proof of Theorem [T] 



□ 



4 The finite dimensional model: proof of Theorem [3] 

We devote this section to describe the proof of Theorem [3| The proofs of the partial results 
stated in this section are deferred to Sections [5H3 

To prove Theorem [3] we need to analyze certain orbits of system (18) given by Hamiltonian 
h(b) in (19). Moreover, there is another conserved quantity: the mass 



M{b) = Y,\b,\'. 

We obtain orbits given in Theorem [s] on the manifold Ai{b) = 1. 

It can be easily seen that on M{b) = 1 there are periodic orbits Tj given by 



bj{t) 



bk{t) =0 for k^j, 



(31) 



(32) 



which in normally directions have mixed type: hyperbolic in some directions and elliptic others. 
Moreover, there exist two families of heteroclinic orbits, which connect consecutive periodic or- 
bits. Consider the 2-dimensional complex plane Lj = {\/k ^ j,j + 1 : 6^ = 0}. In Section [2. 1| we 
show that they are invariant and dynamics inside are integrable. Then, the (two dimensional) 
unstable manifold of the periodic orbit {bj{t), bj+i{t)) = (e"**, 0) coincides with the (two dimen- 
sional) stable manifold of (6j (t), 6j+i(t)) = (0,e^**) and it is foliated by heteroclinic orbits. As 
usual, the stable and unstable invariant manifolds have two branches and, therefore, we have two 
families of heteroclinic connections. It turns out that they can be explicitly computed jCKS"'"10] 
and are given by 



7f(t) = (0,...,0,6,(t),6f+,(t),0,...,0) 



(33) 
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with 

bj{t) = , 6± (t) = ± ^ ^9 G T. 

To prove Theorem [3] we look for an orbit which shadows the sequence of separatrices, as follows 

• it starts close to the periodic orbit T3 

• later it passes close to the periodic orbit T4 

• later it passes close to the periodic orbit T5 and so on 

• finally it arrives to a neighborhood the periodic orbit TAr_i. 
Our main goal is to prove 

existence of such orbits and estimate the transition time in terms of N. 
In making these transition we have the freedom of whether to travel close to 7^ or 7" . We 
will choose always 7^ The procedure for 7" is analogous. 

We believe it is helpful to the reader to have the following information about transition of 



energy. We have a solution b{t) = {6j(t)}j=o,...,Ar to the system (18). We fix a > small, but 
independent of N, and 6 = e~"'^ . For each j = 2, . . . ,N — 1 near the periodic orbit Tj and later 
near Tj^_i we have the following table of orders of magnitude of distribution of energy 

near Tj — > near T^+i 

\b<j-2\ |6<,_2|(l + 0(5'-')) 

\bj-2\ — > J<^\bj-2\ 

\bj.i\=0{a) (C(j)5)i/2 (34) 

\bj\ = 1 — O(cr^) (mass conservation) — > 0(a) 

= {C^^^6)^^'^ — > 1 — 0((T^) (mass conservation) 

\bj+2\ > ^l&j+2| 

\b>j+2\ |6>,+2| (l + 0(5''')). 

We decompose a diffusing orbit into A^ — 5 parts: near each periodic orbit Tj, j = 3, . . . , A''— 1 
we construct sections transversal to the flow so that they divide the orbit appropriately. For 
each transition from one section to the next one we associate a map which sends points close 
to Tj to points close to Tj+i. This leads to analysis of the composition of all these maps 

B* = B^-^ o...oB^. 

To study these maps we will consider different systems of coordinates which, on one hand, will 



take advantage of the fact that mass (31) is a conserved quantity, and on the other hand, will 



be adapted to the linear normal behavior of the periodic orbits. These systems of coordinates 



are specified in Section 4.1 
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4.1 Symplectic reduction and diagonalization 

To study the different transition maps we use a system of coordinates defined in [CKS^IO] . It 
consists of two steps: 



• A symplectic reduction uses that mass (31) is conserved and sends the periodic orbit Tj 
into a critical point. 

• A hnear transformation diagonahzes the hnearization of dynamics near this critical point. 

We perform the change corresponding to the traveling close to the j periodic orbit Tj. We 
restrict ourselves to A4{b) = 1 and we take 

= r(j)e^^''\ 6fc = c[f\*^*'^ for ah k ^ j, (35) 

where 0^^^ is a variable on Tj. From now on in this section we omit the superscripts (j). It can 
be seen that after eliminating r using that Ai{b) = 1 and omitting the equation for the variable 
9, one obtain a new set of equations whose Ck components form a Hamiltonian system with the 
Hamiltonian 

k^j \ k^i ) k^j,j+l 

and the symplectic form Q = ^dck A dck. The Hamiltonian H^^\c) can be written as 

F(^) (c) = H^^'^ (c) + i^l^'^ (c) (36) 



with 



k^j 

k^j \k^j J k^jj+l 

Since we are omitting the evolution of the variable 9, the periodic orbit Tj has become now a 
critical point for the equation associated to this Hamiltonian, which is defined as c = 0. For the 
same reason, the two families of heteroclinic connections defined in (33), now have become just 
two one dimensional heteroclinic connections. 

The second step is to look for a change of variables which diagonahzes the vector field around 
this critical point. This change only modifies the coordinates (cj„i,Cj+i) and is given by 

c,-i \ ^ / iv'pi + uqi \ ^33) 

Cj+i J \ UJ^P2 + U}q2 J 
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where ui = e^'^*/^ (see |CKS"'"10j ) . Note that this change is conformal and leads to the symplectic 
form 

Q = -dck A dck + dpi A dqi + dp2 A dq2- (39) 

To study the Hamiltonian expressed in the new variables let us introduce some notation. We 
define 

Vj={l<k<N;k^j-l,j,j + l}, (40) 

which is the set of subindexes of the elliptic modes. From now on we will denote by q and p all 
the stable and unstable coordinates q = (91,92) and p = {pi,P2) respectively and by c all the 
elliptic modes, namely with k G Vj. 



Lemma 4.1. The change (38) transforms the Hamiltonian (36) into the Hamiltonian 

{p, q, c) = (p, q, c) + H{ (p, q, c) (41) 

with homogeneous polynomials 



H^^\p,q,c) = - ^ |cfc|^ + \/3(pigi +P292j 



and 
where 



Hi'^ (p, q, c) = H\^>^ (p, q) + H^J> (c) + H^^i {p, q, c) 



k=l k=l k,i=0 

ttll!(c)=jT.M' + l(Y.M'] (42) 



-2„2 



kf^j-l,j,j+l,j+2 



Fiji^(p, I.c)=-—Y, |cfc|2 (gipi + q2P2) (43) 



1/9 \ 2 2 / 2 \ 2 2 



1 z', ,2„_ I ,,„_\2^^— 2 1 A ,2„_ , ,,^_\2„2 

2 



- - >2 + ^^92) Cj+2 - - (w 92 + WP2) 



/or certain constants and I'k, Vki G M. 

Remark 4.2. E'uen though the proof of this lemma is a simple substitution of {p, q) we do need 
specifics of the form of the decomposition into Hamiltonians: 

• is the direct product of two linear saddles {pi,qi), i = 1,2 and N — 2 linear elliptic 

points {cfcjfc, k £Vj. 
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• -ff^yp consists of some only saddle terms. In particular, it does not contain terms pf,qf, i = 

1,2 so {q = 0} and {p = 0} are invariant manifolds of H if we set c = 0. This implies 
that the two heteroclinic orbits which connect the critical point {p, q, c) = (0, 0, 0) to the 
next periodic orbit Tj^i are just defined as 

{pf{t),qt{t),pf{t),qf{t),c^it)) = To, 0,-^^,0,0 

Moreover, Tj+i is now defined as \cj^i\ = 1. Due to (38) it is equivalent to P2 + q2—P2Q2 = 
1. 

• Near p = q = 0, which corresponds to the periodic orbit Tj Hamiltonians H^^^ and H^l^ 
are almost integrable. The only source of non-integrability comes from the second line of 
(4-2) for and from the second and third line of (4-3) for H^j^. 

• Later we select regions with c's being exponentially small in N. As the result coupling 
between hyperbolic variables - {p,q) and elliptic ones c's is exponentially small in N. This 
decoupling at the leading order is crucial for our analysis. 

• Among all the constants which appear in the definition of Hamiltonian ( |41[ ), fo2 7^ 
is the only one which plays a significant role in the proof of Theorem Indeed, the 
corresponding term is resonant and will be the leading term in studying the transition 
close to the saddle. We assume, without loss of generality that vq2 > since the case 
uq2 < can be done analogously. 

Proof. To obtain the explicit form of H^''\ note that H^^\c) in (37) can be rewritten as 

^i'^(^) =\ E i^^i' + i E i^^i' + ^i-i + + s'+i + ^1+1 

1 -2 2 I 2-2 ^ i 2 1-2 I 2 , -2 \2 

~ 2 CfcCfc_i + Cj^cj^_i — - [Cj_i + Cj_i + Cj_^^i -t- Cj_^_l) . 

~ (i) 

Written in this way, the second term in the first row is just a constant times H2 squared. 

Then, the particular form of -f^hyp' -^eii ' ^^'^ ^mL obtained just performing the change 

of coordinates. □ 



Since the symplectic form is given by (39), equations associated to the Hamiltonian (41) are 

hyp ^ '^iJi-^mix 

r(i) _ a S-O') 



Pl = V3pi + ^hyp.pi + ^mix,pi = V3pi + dg,H^J^ + dq.Hl^.^ 



qi = -VSqi + ^hyp,<?i + ^mix,5i = -VSqi - dp.H^J^ - dp.WjJ^ 

P2= ^/3p2 + ^hyp,p2 + ^mix,p2 = V3p2 + dg.Ill^J^ + Og^Ili^l (44) 

q2 = -V3q2 + ^hyp,g2 + ^mix,g2 = -V3q2 - dp^Il'^^^ - dp.H^jJ^ 

Ck = ick + -2^eii,cfc + -2^mix,cfe = ick " 2idc^H'£ - 2idc^^H^}^. 



25 



where 



-2^hyp,pi = X^(4 - k)vkp\ql ^ + t^12PipI + J^iiPiP2q2 + I'lQPiql (45) 



k=l 



Z'hjp^qi = - ^ kVkp'l ^ - 2z^22PlPi - 2l/2lPlP2g2 " 2z^20Pl^2 (46) 



k=l 

J^uqipl - J^iiqiP2q2 - J^wqiqi 



^hyp,p2 = X](4 - k)vkplql ^ + ^2ip\p2 + viipiqip2 + i^oiqfp2 (47) 



k=l 

+ 2l^20Pig2 + 2l/ioPl5lQ'2 + 21/00^192 



-2^hyP,92 = - XI ^^'=^^'2 ^^2 ^ - 2f22P?P2 " 2l^l2Pigi^32 " 2f02g?P2 (48) 



T^2ipiq2 - J^iipiqiq2 - T^oiqlq2 



-eii,cfc = - ^Icfcl^Cfc - M X] Ck + 2ick {cl_i + cl^i) (49) 



-2^mix,gi =a;^(a; V + wg'i)cj_2^ + a;(u;pi + w^gi)cj_2^ + k^l^O'i (50) 

-2^mix,pi = - w(wV + u}qi)c~f - u?{u)pi + w^g'i)cj_2^ ~ |q| V (51) 

2^mix,g2 ='^^('^^P2 + a;g'2)c^^ + uj{ujp2 + io'^q2)cj+2^ + ^ X] |c£pQ'2 (52) 

Zram,P2 = - io{u?p2 + a;g'2)^+2^ " W^(tJp2 + Up'q2)Cj+2^ 2~ ^ |c^l^P2 (53) 

i^/3cfc((7ipi + g2P2) ioTkeVj\{3±2} (54) 
^mix,cj_2 =i\/3cj_2(g'iPi + 92^2) - 2i(a; V + LoqifcJZ^ (55) 

2^mix,Cj+2 =«\/3Cj+2(Q'l?'l + 92^2) - 2i{uj^P2 + Wg2)^CH2- 

4.2 The iterative Theorem 

Now that we have obtained the adapted coordinates for each saddle we are ready to explain the 
strategy to prove Theorem [3j To obtain the orbit given in Theorem [3| we will consider several 
co-dimension one sections {S™}^]^ and transition maps from one section to the next 
one ^^j^-^i- Then, we will detect a class of open sets {Vj}j, Vj C S™, j = 1, . . . , N — 1, which 
have a certain almost product structure (see Definition 4.3) such that Vj+i C (Vj) and none 



of them is empty. Each set Vj is located close to the stable manifold of the periodic orbit Tj. 
Composing all these maps we will be able to find orbits claimed to exist in Theorem [3} 
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We start by defining these maps. Tlie first step is to define certain transversal sections to the 
flow. We use the coordinates adapted to the saddle j, {p^^\ q^^\ c^-'^), which have been introduced 



in Section 4.1, to define these sections. Indeed, in these coordinates, it can be easily seen that 



the heteroclinic connections (33), which connect {p^^\ q^^\ c^^^) = (0, 0, 0) with the previous and 



next saddles are defined by (3^,P2^\ c^^'^) = (0,0,0,0) and {p[^\q[^\p''2^\c^^^) = (0,0,0,0) 
respectively. Thus, we define the map from the section 

= {gp) = a} (56) 

to the section 



Here cr > is a small parameter that will be determined later on. In fact, we do not define the 
map in the whole section but in an open set C S™, which lies close to the heteroclinic 
that connects the saddle j — 1 to the saddle j. Then, we will consider maps 

B^:V,CJ:f^ S^i 

and we will choose the sets recursively in such a way that 

V,+i C B^ (V,) . (57) 

This condition will allow us to compose all the maps B^ . Indeed, the domain of definition of the 
map B^'^^ will intersect the image of the map B^ in an open set. 

The sets Vj will have a product-like structure as is stated in the next definition. Before 
stating it, we introduce some notation. We define the subsets of indices Vj in (40), 



P7 = (^ = i.....;-3} 

r+ = {k=j + 3,..^,N]. 

The first set consists of preceding non- neighbor modes to j — 1, the second — of foreseeing 
non- neighbor modes to j + 1. The modes fc = j it 2 are called adjacent. These modes have a 
stronger interaction with the hyperbolic modes. 

Note that we split the non-neighbor elliptic modes in two sets: the -|- stands for future — 
stands for past. Indeed, along orbits we study future modes will eventually become hyperbolic 
in the future, past have already been hyperbolic. Analogously, we call future adjacent — the 
mode c^'^g ^^i^l past adjacent — 

For a point {p^^\ q^^\ c^^'^) e S^", we define c^'^ = (c[^'\ . . . , c'-^) and cj^ = {c'fl^, • • • , ci^^)- 
We define also the projections TT±{p^^\q^^\c^^'^) = c£^ and vrhyp,+ = {p^^\q^^\cf^). 

Definition 4.3. Fix positive constants r S (0,1), 5 and a and define a multi-parameter set of 
positive constants 

= {^^^^'^ -S, MiJ^, -ii -£^P' <p} • (59) 

Then, we say that a (non-empty) setU C has an Ij -product-like structure if it satisfies the 
following two conditions: 
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CI 



where 



and 



U C 



J] X ... X Dj"^ X A/"/ X D^'+^ X . . . X „j 



N 



V±2 



[\c^A<M'£^J^^-''-y^] fork 



C 



Si) 

-j±2 



< M. 



adj,± 



1/2 



/ (n^^) n'^^') li^^ 

1 \Pl ill 1P2 1 '/2 



-C^^h (ln(l/d) + M^M < pf'^ < -C^^H (ln(l/(5) 



p4 . 



(i) 



hyp 



(60) 



(i) 



C2 



where 



and 



M7 X 



y+2 



r 



.(i) 



< m 



(i)r(l-r)/2 
ell " 



} /■ 



'or k G T',^ 



1/2 



M) Ij) JJ) JJ) 
Pi ) yi 1P2 ' '/2 



C(^)(^fln(l/(5) + 



m 



iJ) 
hyp 



< < -C(J')<5 (ln(l/5) 



(i) 
hyp 



(61) 



(i) 



9Xj{P2,q2,(^,S) = 0, Ifi^- 



(i)i 



i4''i<<(c«>*)"'}. 



The function gx {p2,Q2,(^,^) is a smooth function defined in (94). 



(7) 

Remark 4.4. iVote ^/lai /or i/iis product-like sets the variable p^ is selected negative. This is 
related to the fact that 1/92 > (see Remark 4-2)- The reason of the choice of the sign of p^^ 
wil 



ill be clear in Section\^ In particular, see Remark\K 

The domains Vj of the maps will have Xj-product-like structure as defined in Definition 
4.3| Thus, we need to obtain the multi-parameter sets Ij. They will be defined recursively. 



Recall that, to prove Theorem [3j we want to obtain an orbit which starts close to the periodic 
orbit T3. Thus, the recursively defined multi-parameter sets Ij will start with a set X3. 

Definition 4.5. Fix any constants r,r' G (0, 1) satisfying < r' < 1/2 — 2r, K > and small 
5, a > 0. We say that a collection of multi-parameter sets {Xj}j=3^..,^Ar_i defined in (59) is 
{a, 5, K) -recursive if for j = 3, . . . , N — 1 the constants C^^^ satisfy 

< m^ll'^ < m^£^ 
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and all the other parameters should be strictly positive and are defined recursively as 



r' 



'adj ~ 2 '^'^ 

f (i+1) _ KM^^'^ 

'hyp -^^^■'adj, 



The next Theorem defines recursively the product-like sets Vj, so that condition (57) is 
satisfied. 

Theorem 5 (Iterative Theorem). Fix large 7 > 0, small a > 0, and any constants r,r' € (0, 1) 
satisfying < r' < 1/2 — 2r. Then, if we set 6 = e~'^^ , there exist strictly positive constants K 
and C*^^) independent of N satisfying 

<r^K-(^-2), (62) 



and a multi-parameter set X3 ( as defined in (59)) with the following property: there exists a 
{a, 5, K)-recursive collection of multi-parameter sets collection of multi-parameter sets {2j }j=3^.,,^Ar_i 
and Tj -product-like sets Vj C such that for each j = 3, . . . , N — 1 we have 

Moreover, the time spent to reach the section S^P^-^ can be bounded by 

iTg.! < i^ln(l/5) 
for any (p, q, c) G Vj and any j = 3, . . . , N — 2. 
Note that the condition 

implies 

i^-(i-2)(^(3) < (-^(j+l) < K^+^c^3) 

Namely, at each saddle, the orbits we are studying may lie further from the heteroclinic orbit. 
Nevertheless, by the condition on S from Theorem [3] and (62), these constant does not grow too 
much. Indeed, 

6' < C^-') < 5-', (63) 



where r > can be taken as small as desired. We will use the bound (63) throughout the proof 
of Theorem [H 

Theorem [3] is a straightforward consequence of Theorem [5} In fact, we need more precise 
information than the one stated in Theorem [3j This more precise information will be used in 
the proof of Theorem [4| We state it in the following theorem. Theorem [3] is a straightforward 
consequence of it. 
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Theorem 3— bis Assume that the conditions of Theorem hold. Then, there exists an 
orbit b(t) of equations (18), constants K > and ly > 0, independent of N and 6, and Tq > 
satisfying 

To < KNln{l/6), 

such that 

|63(0)|>l-5" ^ \bN-2{To)\ > I - 

and 

|6,(0)|<<5'^ forj^S \b,iTo)\<6'' forj^N-2 

Moreover, call tj G [0, Tq] the times for which b(tj) G S™, Then, 

tj+i - tj < Kln(l/(5) 
and for any t G [tj, tj+i] and k ^ j — 1, j, j + 1, 

\bk{t)\<S''. 



Proof of Theorem 3-bis. It is enough to take as a initial condition 6" a point in the set V3 C Sg" 
obtained in Theorem [5] Then, thanks to this theorem we know that there exists a time Tq 
satisfying 

TQ^N\n{l/5), 

such that the corresponding orbit satisfies that b{TQ) G Va^-i C T}^_^. Note that in this 
section there are two components of b with size independent of 8. Nevertheless, from the proof 
of Theorem [s] in Section [6] it can be easily seen that if we shift the time interval [0, Tq] to 
[pln(l/(5), pln(l/5) + To], for any p < \/3, there exists u > such that the orbit b{t) satisfies 
the statements given in Theorem 3-bis. □ 

4.3 Structure of the proof of the Iterative Theorem [5] 

To prove Theorem [5] we split it into two inductive lemmas. The first part analyzes the evolution 
of the trajectories close to the saddle j and the second one the travel along the heteroclinic 
orbit. Thus, we study as a composition of two maps. 

We consider an intermediate section transversal to the fiow 



yo 



out 



= a] , (64) 
and then we consider two maps. First the local map 

Bi^:V,C^f^ S-S (65) 
which studies the trajectories locally close to the saddle. Then, we consider a second map, 

e^i„b^W'cS-*^S;'Vi, (66) 



which we call global map, that studies how t 
Then, the map considered in Theorem js 



le trajectories behave close to the heteroclinic orbit, 
'is just B^=Bi,^^oBi 
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Before we go into technicalities we write a table analogous to ( 34 ) of the properties of the 
local and global maps. The local map B'^^^, projected onto hyperbolic variables, has the form 












= a 














k?^|<C(^')5 In^ 



(67) 



The global map S-^^i^^, projected onto hyperbolic variables of the corresponding saddles, has the 
form 



\p?\ 














= a 




= a 


\pr''\ 






<C(^) 6 In] 









(68) 



To compose the two maps we need that the set U\ introduced in (66), has a modified 
product-like structure. To define its properties, we consider the projection 



Jj) Aj) Si) 

P2 ' ^2 ' 



Definition 4.6. Fix constants r G (0, 1), 6 > and a > and define a multi-parameter set of 
positive constants 



Then, we say that a (non-empty) setU C S^"* has a Xj -product-like structure provided it satisfies 
the following two conditions: 



CI 



where 



U dB] y. . . . ^ X A/"-; _ X or " X . . . X 



5f 



V±2 



Si) 



<m'£j^'-^'^I^] for kerf 



c 



.(i) 

-j±2 



1/2 



and 



M; = {{p'^^\q^^\pi^\q?)eR': 



P? 



1/2 



pi'^ = <y 



, -C(^) 5 [hi/ 5) + Mg;,) < < -C(^) 6 (ln(l/5) - M^^ }, 
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C2 



Wx 
where 



-C^^U (ln(l/(5) -m^^'^ 



hyp 



X D^:+2 X 



X C n{U) 



J.- 



'or k G T^j 



Si) 



1/2 



With this definition, we can state the following two lemmas. Combining these two lemmas 
we deduce Theorem [5l 

Lemma 4.7. Fix any natural j with 3 < j < N — 2, constants r,r' G (0, 1) satisfying < r' < 
1/2 — 2r and a > small enough. Take 6 = e""'^ , 7 = 7(0") 1, depending on a, and consider 
a parameter set Xj with Afj^^p > 1 and a Xj -product-like set Vj C S™. Then, f or N big enough, 
there exists: 

• A constant K > independent of N and j but which might depend on a. 

• A parameter set Xj whose constants satisfy 

C^^^ /2 < C^^^ < 2C7(J') 



< < m'^l 



and 



W 
hyp 



K 



ii) _ mU) 

adj,± ~ adj,± 



M. 



(1 + 4cj) 



• A Xj -product-like setlAj for which the map B^^^ satisfies 
Moreover, the time to reach the section S^"* can be bounded as 



(69) 



T„j 



< K\n{l/6). 



The proof of this lemma is the most delicate part in the proof of the Iterative Theorem [5| 
since we are passing close to a hyperbolic fixed point, which implies big deviations. It is split 
in several parts in the forthcoming sections to simplify the exposition. First, in Section [5j we 
set the elliptic modes c to zero, and we study the saddle map associated to the corresponding 
system. We call to this system hyperbolic toy model. It has two degrees of freedom. Then, in 
Section [6] we use the results obtained for the hyperbolic toy model to deal with the full system 



and prove Lemma 4.7 



Now we state the iterative lemma for the global maps B"! 



glob' 



32 



Lemma 4.8. Fix any natural j with 3 < j < N — 2, constants r,r' £ (0, 1) satisfying < r' < 
1/2 — 2r and a > small enough. Take 6 = e~"'^ , 7 = 7(0") S> 1, depending on a, and consider 
a parameter set Ij and a Ij -product-like setUj C S^"*. Then, for N large enough, there exists: 



• A constant K depending on a, hut independent of N and j. 

• A parameter set Tj+i whose constants satisfy 



and 



< m\^X'^ < mg)p 



M 



M 



(i+i) 

cU - 

(i+1) 

oll,+ 
'"ell 



max 



ell, 



m 



ell 



KM, 



W 



hyp 



m 



(i+i) 

adj 



^ell 



max 



A Zj-^i -product-like set Vj+i C for which the map BL^^ satisfies 



(70) 



Moreover, the time spent to reach the section S^^-^ can he hounded as 



T„j 



glob 



< K. 



The proofs of this lemma is postponed to Section [7| 

Now it only remains to deduce from Lemmas |4.7| and |4.8| the Iterative Theorem [5| 
Proof of Theorem^ We choose the multiindex X3 so that we can apply iteratively the Lemmas 



4.7 and 4.8, Indeed, from the recursive formulas in Lemma 4.7 and 14.81 it is clear that it is 



enough to chose a parameter set I3 satisfying 



1 < M 



(3) 



ell,- 



« « « ^S,- « ^ell, 



(3) 



and 



we 



From the choice of the constants in X3 and the recursion formulas in Lemmas 4.7 and 4.8 

(i) ^ — ^ — I. — 

have that > 1 for any j = 3, . . . N — 1. This fact along with conditions (^69^ and ^70\i, allow 



us to apply Lemmas 4.7 and 4.8 iteratively so that we obtain the {6, a, i^')-recursive collection 
of multi-parameter sets {Ij}j=3^..,^7v-i and the 2j-product-like sets Vj C S™. In particular, 
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note that the recursion formulas stated in Theorem [5] can be easily deduced from the recursion 
formulas given in Lemmas 4.7 and 4.8 and the choice of X3. 



Finally, we bound the time 



IT, 



< 



+ 



glob 



< {K + K)\n{l/5). 



This completes the proof of Theorem [5] 



□ 



5 The hyperbolic toy model 



In this section we set the elliptic modes to zero, namely, we deal with the system 

Pl = V^pi + 2hyp,pi 

^1 = -V^qi + -2:hyp,gi 

P2 = V3P2 + 2^hyp,p2 

q2 = —V3q2 + -Zhyp,g2, 



(71) 



where the functions 2^hyp,* are defined in H5l mh, M7| and M8L 
We start by setting some notation. We call 

z = {xi,yi,X2,y2) 

the new set of coordinates, whose components are also denoted by Zi = {xi,yi). We also use the 
notation x = (xi,X2) and y = (^1,2/2)- 

Moreover, we call K to any positive constant independent of 5, N, j, and a and we call 
Ko- to any positive constant depending on a, but independent of 5, N and j Analogously, we 
say that a = 0(b) if \a\ < K\b\ and that a = Oaib) if \a\ < Ku\b\. We will also use all these 
notations in Section [6] and Section [7l 

The first step is to perform a resonant normal form in a neighborhood of size a of the 
saddle. Note that we do not need much regularity for the normal form since all our study will 
be done in the norm. It turns out it is enough to consider a normal form. Before we state 
our next claim about the normal form we formulate a well known result of Bronstein-Kopanskii 
|BK92j about finitely smooth normal forms of vector fields near a critical point. We are unable 
to use classical results about linearizability, because our saddle is resonant. 

The main result of Bronstein-Kopanskii |BK92| is that near a saddle point a vector field 
can be transformed into a polynomial one by a finitely smooth change of coordinates with only 
certain (resonant) monomials present. For convenience of the reader we use notations of this 
paper. 



5.1 Finitely smooth polynomial normal forms of vector fields in near a saddle 
point 

Let X = F{x) be a vector with the origin being a critical point, i.e. -F(O) =0, x G M"* for some 
d G Assume that F is for some positive integer K G Z_|_, i.e. F has all partial derivatives 
of order up to K uniformly bounded. Denote the linearization of -F at by A := DF{0) and 
f{x) = F{x) — A{x). Then, the equation becomes 

x = Ax + f{x), /(O) = 0, I)/(0) = 0. 
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Let vi, ... ,1^(1 denote the eigenvalues of A and 9i, . . . ,9n be all distinct numbers contained in 
the set {Rei'i : i = 1, . . . , d}. Assume that none of ^j's is zero or, in other words, the rest point 
being hyperbolic. 

The space M'* can be represented as a direct sum of j4-invariant subspaces i?i , . . . , En such 
that the eigenvalues of the operator AIe^ satisfy the condition Refj = 9i. 

Theorem 6. \BK9^ Let k be positive integer. Assume that the vector field x = F{x) is of class 
, X = is a hyperbolic saddle point and A = DF{0). If K > Q{k) for some computable 
function Q{-), then, for some positive integer N , this vector field near the point x = can be 
reduced by a transformation y = <I>(x),<I> G , to the polynomial resonant normal form 



where r G Z'j^ andpr denotes a multi-homogeneous polynomial Pr{Ei, . . . , EiQ)- ■ -(BEn), Pt = 
{p]j., . . . ,p'^) and / implies Vi = t^vi + • • • + r'^Vd (by the resonant condition). 

In Theorem 3 |BK92j the authors give an upper bound on N . In our case d = 4, n = 2, k = 1. 
A direct application of this Theorem is the following 

Lemma 5.1. There exists a change of coordinates 



N 




\t\=2 



{Pl,qi,P2,q2) = '^hyp{xi,yi,X2,y2) = (xi , ?/l , X2 , ^2) + ^'hyp(xi , 7/1 , X2, 2/2) 



which transforms the vector field ( 71 ) into the vector field 




where D is the diagonal matrix D = diag(\/3, — \/3, V^i — \/3) and Rhyp is a polynomial, which 
only contains resonant monomials^ It can be split as 




and -R^vn ^■^ remainder and satisfies 





O {yl,xiyi,yi{xl + yl), xiX2{x2 + ^2)) 




^hyp,j/2(2;) = O {yl,x2y2,y2{x\ + yl), xiX2{xi + yi)) . 
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One can even estimate degree of this polynomial using |BK92] 
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5.2 The local map for the hyperbolic toy model in the normal form variables 

Recall that our goal in this step of the proof is to study the evolution of points with initial 
conditions inside of a certain set near the section S™. More specifically, in formulas (60) and 
(61 ) we define sets Mj' C A^^. We set elliptic modes c = and shall study the set AAj satisfying 

AAr n {c = 0} C A/} C AA+ n {c = 0}. 

Since the analysis is done in normal coordinates ^hyp : {x,y) — ?• {p,q), we study the a set Afj 
such that ^hyp(-^j) ^ -^r define this set we need to fix several parameters and define several 
objects. 



Let C'--'-' 's be the constant from Lemma 4.7 Recall that in Definition 4.5 we define a (a, (5, K) 



recursive multiparameter set Ij. Its description includes parameters Af^yp used below. The 
parameter K depends on a and we keep this dependence in the notation: K^j. Denote the 
inverse of the map ^' from Lemma |5.1[ by 



T := Id + T := ^^^^ =: Id + (T,, , T,.„ T,,). 

Define 

a(^) := C(^) (l + a,,f,,(0,c7,0,0)) . (76) 

Notice that d^^^ = C^^^ (1 + 0{a)) . Define fi{(T) by 

/i((t) = Tj,, (0,(7, 0,0). (77) 

Observe that it satisfies /i(o") = a + 0{a^) and the section {yi = /i(cr)} approximates the image 
of the section T(S™). Now we can define the set of points whose evolution under the local map 
we shall analyze 



Mj = {|xi + C(^')5(ln(l/(5)| < d^^UK„, \x2-x*2\ < 2Mj 



(i) 



ln(l/5) ' (78) 
\yi - fi{<T)\ < K.a(^)<^ln(l/<5), \y2\ < 2MgJ, (a(^)<5)'^' }, 

where the constant X2 will be defined later in this section. It turns out a proper choice of X2 
leads to a cancelation in the evolution of the xi coordinate (described in Section 2.2 for the 
simplified model). This cancelation is crucial to obtain good estimates for the map B'^^^. 
We also define the function /2(o") as 

f2{<7) = T,,{0,0,a,0). (79) 

By analogy with /i(o") notice that the section {x2 = /2(o')} approximates the image of the 
section T(S°'^*) with 5]°^* = {p2 = a}- Later we need to compute an approximate transition 
time Tj{x2) from near T(SJ'^) to T(S^"*). We use /2 to do that. Notice that the X2 coordinate 
behaves almost linearly as 

X2 ~ x°e^*. 

Therefore, for an orbit to reach {x2 = /2(o")} it takes an approximate time 
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Note that this time is defined for any X2 > 0. We will see that the X2 coordinate behaves as 
2 — (C^^')6y^'^ and, therefore, Tj behaves as 

1 



In 



C(i)5 

Even if X2 behaves approximately as for a linear system, this is not the case for the other 
variables, as we have explained in Section 2.2 with a simplified model. Indeed, if one first 
considers the linear part of the vector field (71), omiting the dependence on C^^\ the transition 
map sends points 

(xi, yi, X2, y2) ~ [0(5 Hl/6)), Oia), O {s^'^) , O {s^'^) ) 



(xi, yi, X2, ys) ~ [O (d^l^ Hi/ 5)) , O (d^'^) , 0{a), 0{5) 



to 



However, the resonance implies a certain deviation from the heteroclinic orbits. Indeed, one can 
see that tipically, the image point is of the form 

(xi,yi,X2,y2) ~ ip {5^'Hn{l/5)) ,0 {^5^'^) ,0{cj),0{6Hl/5)) . 

This apparently small deviation, after undoing the normal form, would imply a considerably big 
deviation from the heteroclinic orbit and would lead to very bad estimates. Nevertheless, if one 
chooses carefully X2 in terms of xi and yi, one can obtain a cancelation that leads to an image 
point of the form 

(:ri,yi,X2,y2) ~ (O (5'/^),0(5'/^),0ia),0{6\n{l/5) 



Since the points we are dealing with belong to the set Mj defined in ( 78 ) , this cancellation boils 
down to choosing a suitable constant X2- Next lemma shows that a particular choice of X2 
leads to a cancellation that allow us to obtain good estimates for the saddle map in spite of the 
resonance. The choice we do is essentially the same as the one choosen in Section 2.2 for the 
simplified model that has been considered in that section. 

Lemma 5.2. Let us consider the flow ^^-^^ associated to (72) and a point G J\fj. Then, if we 
choose X2 as the unique positive solution of 



{xlYT.ixD 



C(J),5 ln(l/(5) 
2 1^02/1(0-) 



(81) 



and we take 6 and a small enough, the point 



where Tj = Tj{x2) is the time defined in (80), satisfies 



I /i 

mi 



14-/2(^)1 
y/ + AMao)5in(i/5) 



1/2 
1/2 
1/2 



H{l/5) 
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Remark 5.3. The particular choice of X2 being a solution (81) will ensure a cancellation. This 
cancellation is crucial to obtain good estimates for the local map. 

Equation (81) has real solutions becau se vq2 > (see Remark 4-2) and xi < (and pi < 
in the original variables, see Remark 4-4^- Indeed, if xi > and xi ~ 0(^)5 ln( 1/(5) we have 

(x*?T (x*)- ^^'^ ^ ^"^^/^^ 

If there is no solution to this equation, we cannot attain the desired cancellation. 

Let us point out that taking into account the estimates for the points in Af^^\ the definition 
of Tj in (80) and condition (63), one can deduce that condition (81) imphes 

1/2 



and then, 



Tj{x°2)<K^\n{l/5). 



(82) 

I I f f 

We use this estimate throughout the proof of Lemma 5.2 Note also that for the modes {x-^ , y( ) 
we just need upper bounds, since after the passage of the saddle j, the associated mode will 
become elliptic and therefore we will not need accurate estimates anymore. 



Proof of Lemma 5.2. We prove the lemma using a fixed point argument. We look for a con- 
tractive operator using the variation of constants formula. Namely, we perform the change of 
coordinates 

Xi = e'^'ui, Vi = e-'^'vi (83) 
and then we obtain the integral equations 



Ui 



dt 



(84) 



dt. 



In the linear case Uj's and Uj's are fixed. We use these variables to find a fixed point argument. 



We define the contractive operator in two steps. This approach is inspired by Shilnikov [Sil67 
First we define an auxiliary (non-contractive) operator we follows 



.^hyp — (.^hyp,tii ; ■I'hyp,Vi ; •^hyp,M2 ! -^hyp,' 



V2J 



as 



•^hyp,n,(w,w) = + /" 

Jo 



ue^Kve-^'' 



dt 
dt. 



(85) 



One can easily see that in the ui and V2 components the main terms are not given by the initial 
condition but by the integral terms. This indicates that the dynamics near the saddle is not well 
approximated by the linearized dynamics and the operator is not contractive. 
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We modify slightly two of the components of J^hjp and obtain a contractive operator. We 
define a new operator 

•^hyp — (-^hyp,Jii ) -^hyp,?)! 1 •^hyp,U2 ' '^hyp,t)2 ) 



as 



(86) 



-^hyp,«i iui,Vl,U2,V2) = ^hyp,ui {ui , -7^hyp,i;i ("Ul , t^l , "2 , ^'2 ) , -^hyp,«2 (^1 , , ""2 , t^2) , ^2 ) 

-^hyp,i;i(Ml,fl,n2,f2) = -7liyp,t,i (u-l , , ""2, t'2) 

-^hyp,«2("l;^^l;'"2,W2) = -7liyp,u2 (^1 ' ^1 ' ^2 , 1^2) 

-^hyp,i;2 (^1 , ""1 , ""2 , ^2) = -^hyp,D2 (^1 > -^hyp.i-i (^1 , t^l , n2 , ^2) , -^hyp.na (""l , 1 , ^^2 , ^'2) , ^^2) 

Note that the fixed points of these operators are exactly the same as the fixed points of -Tiiyp- 
Thus, the fixed points of the operator Jiiyp are solutions of equation (84). 

It turns out the operator J^hyp is contractive in a suitable Banach space. We define the 



following weighted norms. To fix notation, we denote by 
Then define 



the standard supremum norm. 



hyp, Ml 



\h 



hyp, 1)1 
hyp,M2 

^11 hyp, 1)2 



sup 

te[o,rj 



C(^)51n(l/(5) +21/02/1 ((t) [xlft + C^^^S] h{t) 



(87) 



: (X 



\h\\ 



and the norm 



[u,v) 



sup i||Ui||hyp,«, 

■t=l,2 



Ihyp,?;;} • 



(88) 



This gives rise to the following Banach space 

3^hyp = {(^x,^;):[0,r] 



[u,v\ 



I* < 00} 



The contractivity of -Tiiyp is a consequence of the following two auxiliary propositions. 



Proposition 5.4. Assume (81), then there exists a constant > independent of a, 6 and j 



such that for 6 and a small enough, the operator J'hyp satisfies 



Proposition 5.5. Consider w,w' G B{2kq) C J'hyp and let us assume (81), then taking 6 <^ a, 
the operator J^hyp satisfies 

1 /2 

WAypiw) - Ayp{w')\U < K„ ip^''^^) \n\l/5)\\w - w'\U. 

These two propositions show that -T^iyp is contractive from B(2kq) C 3^hyp to itself. Moreover, 
using them we can deduce accurate estimates for the image point. We prove here Proposition 



5.4 The proof of Proposition 5.5 is deferred to the end of the section. 
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Proof of Proposition 5.4 ■ We bound each mode separately. For J-hyp,^^ and J-hyp,n2) we have 
that 

-^hyp,«i (0) = Vi and -7iiyp,„2 (0) = 

and therefore, they satisfy the desired bounds. Now we bound the first iteration for ui. Here we 
use the particular choice of in terms of (x^, y?) done in (81 ) to obtain the desired cancellations 
(see Remark 5.3). Indeed, taking into account the properties of -Rhyp,xi given in Lemma 5.1, the 
first iteration is just 

-^hyp,ni(0)(t) =X? + I' {2uo2yl{xlf + OHy'ifixlf) dt 



Therefore, taking into account that G Mj (see ([78])) and also (82), we have that 
Ayp,^.i(0)(t) = + 2uo2fi{cj){xl)H + O ic^^h] 



Thus, applying the norm given in (|87|), we have that there exists a constant kq > such that 



-^hyp,Mi (0) 



hyp, Ml 



To bound the first iteration for V2., we just have to take into account that it is given by 

ft 







2u,2xl{y'lY + O [{y',f {xlf 



^hyp,.2(0)(t) = y2° 
Then, recalling that G AAj, 

-^hyp,.2(0)(t) <^vo2xl{y\fTj, 



dt. 



which gives 

Therefore, we can conclude that 



^(0) 



hyp, 1)2 



< 4f02- 



< 



for certain constant > independent of 5, a and j. 



□ 



The previous two Propositions show that ^hyp is contractive from B{2kq) C 3^hyp to itself. 
Therefore, it has a unique fixed point in B{2kq) C ^ hvp which we denote by w* . Now it only 
remains to deduce the bounds for z^ stated in Lemma 5.2 To this end, we use the contractivity 
of the operator J-hyp and we undo the change ( 83 ) . Using the definition of Tj in ( 80 ) , we obtain 



f^^^> + J-hyp,.,(«;*)(T,) - J-hyp,.,(0)(T,) 



X. 







--f2{a) (^l + o(^(aa(^)(5)'^'ln^ 
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Analgously, one can see that 



1^1 



1/2 



f 

To obtain the estimates for x\, note that the particular choice that we have done for X2 in (81) 
imphes that 

\ui{Tj)\ < |-Fhyp,«i(0)(T,-)| + |Ayp,«iK)(T,-) - Ayp,«i(0)(T,-) 



Then, undoing the change of coordinates (83) and using the definition of Tj in (80), one obtains 



Finally, proceeding analogously, and taking into account (81) again, one can see that 



,ln(l/i) 



which completes the proof of Proposition 5.2 



□ 



Now, it only remains to prove Proposition 5.5 



Proof of Proposition 5.5. To compute the Lipschitz constant we need first upper bounds for 
w € B{2kq) C 3^hyp in the classical supremmum norm || * ||oo* They can be deduced from the 
definition of the norms || • ||hyp,* in (87) and the fact that z'^ G M'^^'^ (see (78)). Then, we have 
that 



\ui\ < K^C'^^H\n{l/5) 

< KCT 

1/2 



\U2 



(89) 



\V2\ < K„ (^C(j)(5)^^^n(l/(5). 



where i^T > is a constant independent of cr. 

We use these bounds to obtain the Lipschitz constant. We start by computing the Lipschitz 
constant of Jiiyp,i;i = -^hyp.i;! and Jiiyp,?i2 = ^'iiyp,u2 ^ud then we will compute the other two. 

Using the properties of i2hyp,2/i given in Lemma 5.1, (82) and the just obtained bounds, one 
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can easily see that 

|-7iiyp,t,i(u,u) - -Fhyp,t,i(n',f')| < / 0{uv) \vi-v[\dt+ I O [v^) \ui - u[\dt 

•^0 i=l,2 -^0 i=l,2 



<K^ {c^^hy\n{l/5) ^ I.e., - .,1100 

i=l,2 



\Vi - V," 



+ K,ln(l/,5) ^ 



\Ui — u. 



illoo 



i=l,2 



1 /2 

<i^^ (c^-'^t^) ln(l/5) II 

1=1,2 

1 /2 

+ (0(^)5) ln(l/5) ^ II 

''^i ''^i II hyp, Mi ■ 

1=1,2 

Note that we are abusing notation since inside the 0{-) the dependence of the size on {u,v) 
means both dependence on {u,v) and {u',v'). We do not write the full dependence since both 
terms have the same size. Applying the norms defined in (87), we get 

<KjC^^)5) ln(l/5)||(n, 
Now we bound the Lipschitz constant of J^hyp,M2- Proceeding as in the previous case one obtains 

\Thyp,U2{u,'v) - ^hyp,U2{u ,v')\ < / O (uv) \Ui - Ui\dt + / O {u^) \Vi - Vi\dt 

■^0 *=1,2 i=l,2 



i=l,2 

Vjlloo- 



\Vi - v'A 



I _ 'II 



j=l,2 
1=1,2 

+ K„d^^^5Hl/6) Y ll^i-^'^llhyp,., 
i=l,2 

and thus 

n1/2 ^ ^ 

||-^hyp,M2 ^) ~ -^hyp,n2(''^ )|lhyp,M2 

<K„lC^i)5] ln(l/5)||(^x , v) — {u\ v') 



To bound the Lipschitz constant of -Fhyp,«i we use its definition in (86). First we study 
F\xyp,ui{w) — -^hyp,Mi (^')- We proceed as for -Fhyp,u2 but we have to be more accurate. We 
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obtain 



f-Tj j-Tj 
|-^hyp,«i('",^') --^hyp,«i I < / 0{UV) \Ui-u[\dt+ I O [u^) 

•^0 i=l,2 i=l,2 



Vi — v[\dt 



1/2 _ 



i 



Ui - u. 



i=l,2 



j=l,2 



1 /2 

<i^. (d^^^5) d^^^5ln\l/6)\\ui - ui\\tjp,u, 
+ K^d^^^5Hl/6)\\vi - v[Ujp,,, 

1/2 

+ K„ (a(^')<5) C(^^6ln\l/6)\\V2 - ^^||hyp,.2- 
Thus, taking into account that for 6 smah enough, 



sup 

te[0,T,(xO)] 



-C7(i)(51n(l/5) + 2z^o2/i(f^) (4) i + C(j)(5 



< 



CO-) 5' 



one can deduce that 



1 /2 

||-7liyp,«i(^i,^^) - -^hyp,«iK,^^')||hyp,«i (^C^^^S^ ln^{l/6)\\ui - ni||hyp,„i 

+ K^ln{l/6)\\u2 - n'allhyp.ua 
+ K„ln{l/6)\\vi - v[\\hyp,v^ 

1 /2 

+ (a(^')<5) ln2(l/5)||t;2 - f^||hyp,.2- 

Therefore, to obtain the Lipschitz constant for J-hyp,un it only remains to use its definition in 
(86) and the Lipschitz constants already obtained for Jiiyp,s;i and J-'hyp,u2 to obtain 

1/2 

Ayp,uAu,v)-Ayp,uAu',v') < K„ (c^^^d) ln\l/6mu,v) - {u',v')\U. 

hyp,«i V / 

Proceeding analogously, one can see also that 

-^hyp,.2K^') --^hyp,.2(n',^;') < (C^^hy^\n{l/5mu,v) - {u',v')\U. 

hyp,ti2 V / 

This completes the proof. □ 



6 The local map: proof of Lemma 4.7 



Analysis of Section [5| describes dynamics of the hyperbolic toy model ( 71 ) . Now we add the 
elliptic modes and consider the whole vector field (44). Our goal is to study the map B^^^. The 



key point of this study is that the elliptic modes remain almost constant through the saddle 
map and do not make much influence on the hyperbolic ones. In other words, there is an almost 
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product structure. This allows us to extend the results obtained for the hyperbolic toy model 
(71) in Section [5] to the general system. 

As a first step we perform the change obtained in Lemma 5.1 by means of a normal form 
procedure for the hyperbolic toy model (71). The proof of this lemma is straightforward taking 



into account the form of the vector field (44) and the properties of ^'hyp given in Lemma 5.1 



Lemma 6.1. Let ^'hyp be the map defined in Lemma 5.1. Then an application of the change of 
coordinates 

{Pi,qi,P2,q2,c) = {'ifhyp{xi,yi,X2,y2),c) , (90) 



to the vector field (44) leads to a vector field of the form 

Z = Dz + i?hyp(^;) + Riai^,z{Z: c) 
Cfc = iCk + 2:ell,cfe(c) + i?mix,c(2:, c), 

where z denotes z = {xi,yi,X2,y2), D = diag(-v/3, — -v/3, -v/3, — \/3), -^hyp has been given in 

are defined as 



Lemma 



5.1 



Z'c\\,ck is defined in (49), and -Rmix,^ and Rmix,ck 



R 



miXjXi 



A^^{z)cJZ^'^ + A^^{z)cj.2'^ + ^Y1 



R. 



■mix, 1^1 



Ay^{z)cjzf + Ay^{z)Cj_2^ + ^ 



R 



mix,X2 



^^•2(^)9+2^ + A:,^{z)Cj+2'^ + ^ X] |Cfc|^^X2(^ 



kev 



R. 



mix, 1^2 



Ay^{z)Cj + 2 



kep 



R 



miXjCj. 



R 



mix,Cj±2 



i^/3ckP{z) for m ^ j ±2 
iV3cj±2P{z) - icj±2Q±{z) 



where ^hyp,z are the functions defined in Lemma 5.1, A^ satisfy 

Axi = 0{xi,yi) and Ay- = 0{xi,yi) 

and P and Q± satisfy 

P{z) = O {xiyi,X2y2,zlzl) , Q^{z) = O {xi,yi) and Q+{z) = O {x2,y2) ■ 
One can easily see that for this system there is a rather strong interaction between the 



hyperbolic and the elliptic modes due to the terms -Rn 



and R 



mix,yi • 



The importance of 



these terms can be seen as follows. The manifold {x = 0, y = 0} is normally hyperbolic |Fen741 
IFen771 IHPS77j for the linear truncation of the vector field obtained in Lemma 6.1 and its 
stable and unstable manifolds are defined as {x = 0} and {y = 0}. For the full vector field, 
the manifold {x = 0,y = 0} is persistent. Moreover it is still normally hyperbolic thanks to 
|Fen74l IFen771 IHPS77] . Nevertheless, the associated invariant manifolds deviate from {x = 0} 
and {y = 0} due to the terms Rmix,Xi and Rmix,yi- To overcome this problem, we slightly modify 
the change (|90[) to straighten these invariant manifolds completely. 
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Lemma 6.2. There exist a change of coordinates of the form 

iPi,qi,P2,q2,c) = (^(xi,yi,X2,y2,c),c) = (xi,yi,X2,?/2,c) + (^$(xi,yi,X2,?/2,c),0^ 



which transforms the vector field ( 44 ) into a vector field of the form 

Z = Dz + -Rhyp(z) + -Rmix,^(^, c) 
Cfc = iCk + ^eW^ck (c) + Rma,Ck {z, c) , 



where Rhyp cii^d Z^w are the functions defined in (73) and ( |49| ) respectively, and 



^mbc,2/i = 52;i(2,c)c^^ +Sj^i(z,c)Cj_2^ + ^ X] I ^fc I ^C"?/! (2, c) 

fceP 

-Rmix.cfe = iVSckPiz, c) forkj^j±2 

^mix,Cj±2 = i\/3Cj±2^'(^, c) - «Cj±2Q±(^,c), 

where the functions and Cz satisfy 

B^^ (z, c) = (xi + yi2;2^;2) B^^ {z, c) = O (x2 + ^23:1^1) 

By^ {z, c) = (yi + xiy2Z2) By^ (z, c) = O {y2 + X2yizi) 

(z, c) = {xi + yi 2:2^2) {z, c) = {x2 + y2Xizi) 

Cy^{z,c) = 0{yi + xiy2Z2) Cy^{z,c) = 0(2/2 + X2yizi) 

and P and Q± satisfy 

P{z,c) = O {xiyi,X2y2,zlzl) , Q^{z,c) = O {xi,yi) and Q+{z) = O {x2,y2) 

Moreover, the function ^' satisfies 

^xi=C> i xl,xiyi,xi{xl + yl),yiy2{x2 + y2),c^_2yi, ^ \ck\'^yiyl 1 
\ kev / 

$1/1 = O I yl,xiyi,yi{xl + yi),xiX2(x2 + y2),cj_2Xi, ^ |cfcpxix^ | 
^x2 = o{ xl,X2y2,X2{xl + yl),yiy2{xi + yi), c^+22/i> ^ 



= ( yLx2y2,y2{x\ + yf),xiX2(xi + yi), 0^+22^1. X] I'^'^'l 
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2 2 
X2X1 



Proof. It is enough to compose two change of coordinates. The first change is the change (91) 
considered in Lemma 6.1 The second one is the one which straightens the invariant manifolds of 
a normahy hyperbolic invariant manifold [ Fen74| FFenT?! lHPS77j . Then, to obtain the required 
estimates, it suffices to combine Lemmas 5.1 and 6.1 with the standard results about normally 
hyperbolic invariant manifolds. □ 

After performing this change of coordinates, the stable and unstable invariant manifolds of 
{x = 0,y = 0} are straightened. This will facilitate the study of the transition map close to the 
saddle. 



As we have done in Section |5| we define a set Vj such that 



(93) 



where Vj is the set defined in Lemma 4.7 and T is the inverse of the coordinate change ^ given 
Then, we will apply the flow <I>* associated to the vector field (92) to points in 



6.2 



'J 

in Lemma 

Vj. To obtain the inclusion (93) we define the function gi {p2,Q2,cr,5) involved in the definition 
Define the set 



ofV,. 



X . . . X ^ xMj X D^:+^ X ... X 



^3 ' 



where Mj is the set defined in ( 78 ) and D^' are defined as 



O,^ = { Icfel < Meii,±5(i-^')/2} for k G Vf 
Df2 = ||c,±2|<Madj,± [d'-^hf^Y 

Define the function gXj{P2^(l2^(^-,^) involved in the definition of the set Vj as 

gXj{P2,q2,<y, ^) = P2 + ap{a)p2 + ag{a)q2 - 



(94) 



where X2 is the constant defined in (81) and 

aqicr) = 9g2Tp2 (0,0-, 0,0,0), 
where T = Id + T is the inverse of the change ^ given in Lemma 



6.2 



Lemma 6.3. With the above notations for 5 small enough condition (93) is satisfied. 



Proof. It is a straightforward consequence of Lemmas |5.1| and 6.2 



□ 



After straightening the invariant manifold, next lemma studies the saddle map in the trans- 
formed variables for points belonging to Vj. 



Lemma 6.4. Let us consider the flow associated to (92) and a point (z^, c ) G Vj. Then for 
5 and a small enough, the point 



z-^ , c^ 



46 



where Tj = Tj{x2) is the time defined in {80], satisfies 



1^1 1 



1/2 



1/2 



J 



/2(^)l 



< AM 



6. 



and 



i - 

S±2 S±2^ 



<i^^5(i-^)/2+'^' for k G Vf 



1/2 



We postpone the proof of this lemma to Section 6.1 



Now, to complete the proof of Lemma 4.7 we need two steps. 



The first is to undo the change of coordinates performed in Lemma 6.2 to express the 
estimates of the saddle map in the original variables. 

The second step is to adjust the time so that the image belongs to the section T,"^^. These 
two final steps are done in the next two following lemmas. 

Concerning the first step, recall that the change of variables ^ defined in Lemma 6.2 does 
not change the elliptic variables, and therefore it only affects the hyperbolic ones. 



Lemma 6.5. Let us consider the flow associated to (44) and a point {p ,q ,c ) G Vj. Then 
for 6 and a small enough, the point 



p^ , 



where Tj is the time defined in (80), satisfies 



\P{\ 



1/2 



1/2 



\P2 - ^1 



< Ky 

+ C(j')(51n(l/5)| < C^^^ 6 K^. 
for certain constant C^^^ satisfying C^^^^ /2 < C^^^ < 2C^^'^ and 

1/2 



J ^0 iTj 

S±2 ~ S±2^ 



Proof. In Lemma 6.2 we have defined the change ^' which relates the two sets of coordinates by 



f f f f f f 

[Pl,qi,P2^Q2^C-' 



{x{,y{,x{,y{,i 
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Then, taking into account the properties of the change ^ stated in this lemma, one can easily 
see that from the estimates obtained in Lemma 16.41 one can deduce the estimates stated in 



Lemma 6.5 First recall that the change ^' does not modify the elliptic modes and therefore we 
only need to deal with the hyperbolic ones. 

Using the properties of ^ and modifying slightly , it is easy to see that for 6 small enough. 



\p(\ <K^ (C'-^^S 



m 



{ I <Ka ( c^^h 



1/2 



1/2 



To obtain the estimates for p2 it is enough to recall the definition of /2(o") in (79). For the 
estimates for it is enough to see that from the properties of ^' and the estimates for one 
can deduce that 

52 = (0, 0, CT, 0)X2 + (C^^h^^ 

Therefore, we can define a constant such that the estimate for q2 is satisfied. □ 

Once we have obtained good estimates for the approximate time map in the original variables, 
we adjust it to obtain image points belonging to the section S^*^*. 

Lemma 6.6. Let us consider a point {p^,q^,c^) £ $"^^(Vj), where <l>* is the flow of (44), Tj is 
the time defined in (80) and Vj is the set considered in Theorem^ 

Then, there exists a time T' , which depends on the point {p^,q^,c-^), such that 

{p\q\c*) = ^'^' (pf,qf,cf)&Y.f\ 

Moreover, there exists a constant Kfj such that 

\T'\ < K^5'' 

and 

4 



(95) 



J 



<KJ^-' former 



* f 
Pi - Pi 

* f 
Qi - Qi 



■1— r 



■1— r 



<l2 



P2 

4 



a 



< i^<,C(^)(52-Mn(l/5). 



Proof. The proof of this Lemma follows the same lines as the proof of Proposition |7.3[ Namely, 
first we obtain a priori bounds for each variable, which then allow us to obtain more refined 
estimates. □ 



To finish the proof of Lemma 4.7, we define Uj = Bi^^{Vj) and we check that this set has a 



Ij-product-like structure for a multiindex Ij satisfying the properties stated in Lemma 4.7 (see 



Definition 4.6). Indeed, from the results obtained in Lemmas 6.5 and 6.6 and recalling that by 



the hypotheses of Lemma 4.7 we have that M^^^ > 1, it is easy to see that one can define a 
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and the constant C^^'^ given in Lemma 
in Definition 14. 6[ 



6.5 



constant Kf^ so that if we consider the constants M^^^_^, M^^i j_ and Mj 



r'(i) 



hyp 



defined in Lemma 



4.7 



the set Ui 



^Lci^j) satisfies condition CI stated 



Thus, it only remains to check that the set Uj also satisfies condition C2 of Definition 4.6 



First we check the part of the condition C2 concerning the elliptic modes. Indeed, from the 

one can 

~ (7) 

easily see that for any fixed values for the hyperbolic modes, if one takes the constants rnj^{, 



estimates for the non-neighbor and adjacent elliptic modes given in Lemma 6.5 and 6.6 



m 



(i) 

adj Si'^sn in Lemma 



4.7 



the image of the elliptic modes contains disks as stated in Definition 
4.6 Then, it only remains to check that the inclusion condition is also satisfied for the variable 
52- Prom the proof of Lemma 6.4 given in Section [6.1[ one can easily deduce that the image 



in the 1/2 variable contains an interval of length 0{C^^'6) and whose points are of size sn ialle r 
than 2C^^^ 6 \n{l / 5) . Then, when we undo the normal form change of coordinates (Lemma 6.5) 
this interval is only modified slightly but ke eping still a length of order (C^^^ 6). Thus taking 
into account the constant C^^^ given Lemma 



6.5 



and the results of Lemma 



6.6 



we can obtain a 



constant 'Ti^yp so that condition C2 is satisfied. 

Finally, it only remains to obtain upper bounds for the time spent by the map B^^^. To this 
end it is enough to recall that the time spent is the sum of the time Tj defined in ( |80[ ), which 
has been bounded in (82), and the time T' given in Lemma 



6.6 



which has been bounded in 

(|95|). Thus, taking i nto a ccounts these two bounds we obtain_the bound for the time spent by 

This finishes the proof of Lemma 



^ioc given in Lemma 



4.7 



4.7 



6.1 Proof of Lemma 16.41 

As we have done in the Section[5| we make variation of constants to set up a fixed point argument. 
Namely, we consider 

Xi = e'^^Ui, Vi = e-'^^i, Ck = e'^dk 
and then we obtain the integral equation 



Ui = Xi + I e 

'0 



' (r, f.,^V3t „^-V3t 



''hyp,Xi 



ue^ ,ve 



de" dt 



Vi = y'i + 



yp,yi 



ue^'\ve-^^']+R 



mix,j/i 



dk = 4+ e-'* (Z,n,c, {de'') + R 



de" dt. 



ue^Kve-^\de'']] dt 



(96) 



Note that the terms -Rhyp,z are the ones considered in Section [5| and, therefore, we will use the 
properties of these functions obtained in that section. We use the same integration time Tj in 



(80). 



As before, we use (96) to set up a fixed point argument in two steps. First we define 
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G = (^hyp,^ell) as 

Qhyp,vSu,v,d) = yO - ' e^* (iihyp.y, (^xe^*,7;e-^*) + i?,nix,., (^xe^*, 7;e-^*, de^*)) dt 
= J'hjp,vAu,v) + ^ ' e^*i?mix,x, (ue^',ve-^\de''^ dt, 
where J-^hyp is the operator defined in ( 85 ) , and 

gen,c,{u,v,d) = 4 + J^' e-^* (^eii,c, (de^*) + ^mix,c, (txe^*, t;e-^*, de^*) ) dt. 

We modify this operator sHghtly as we have done for J-hyp in Section [5] to make it contractive. 
We define 

^hyp.Mi {U, V, d) = ^hyp,«i {Ul,ghyp,vi {u, V, d), ^hyp,«2 V, d),V2, d) 
Ghyp,V2 {u,V,d) = ^hyp,i>2 (""l ) ^hyp.-Ui {u,V,d), Qhyp,U2 {u, V,d),V2,d). 

We denote the new operator by 

G = ^^hyp,?ii ) ^hyp,n2 ) ^hyp,iii ) ^hyp,D2 ) ^ell^ j i^'^) 

whose fixed points coincide with those of Q. 

We extend the norm defined in (87) to incorporate the eUiptic modes. To this end, we define 

and 

\\{u,V,d)\\^ = sup ] ||nj||hyp,«,, ||fi||hyp,D,, ||dfc||ell,±, IMi±2||adj,± f 
i=l,2 



I oo 

,-1/2 



which, abusing notation, is denoted as the norm in (88). We also define the Banach space 
3^ = {{u,v,d) : [0,T] ^ X 1^4. ||(^^^^^)||^ < . 

Proceeding as in Section [5j we state the two following propositions, from which one can easily 
deduce the contractivity of Q. The p roof of the first one is straightforward taking into account 
the definition of G and Lemma 5.4 and the proof of the second one is deferred to end of the 
section. 
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Proposition 6.7. Let us consider the operator Q defined in (97). Then, the components ofQ{0) 
are given by 



^hyp,ui (0) 


— J~hyp,Ui (0) 


^hyp,t.i (0) 




^hyp,«2 (0) 


_ ™0 
— X2 


^hyp,t.2 (0) 




^ell,c,(0) 





Thus, there exists a constant ki > independent of a, 5 and j such that the operator Q satisfies 



Proposition 6.8. Let us consider wi,W2 € B{2ki) C y, a constant r' satisfying < r' < 
1/2 — 2r and 6 as defined in Theorem^ Then taking a small enough and N big enough such 
that < 5 = e~'^^ <C 1, there exist a constant > which is independent of j and N, but 
might depend on a, and a constant K independent of j, N and a, such that the operator Q 
satisfies 

^hyp,n, (■", V, d) - Ghyp,u, {u , v' , d!) < 

hjp,Ui,Vi 

< \\{u,v, d) — {u ,v' , d')\\^ 

Qhyp,Vi {U, V, d) - Qhyp,v, {u , v' , d') < 

hyp,Ui,Vi 

< \\{u,v, d) — {u ,v' , d')\\^ 

Qell,Ck (u, V, d) - Qe\\,Ck {u , V , d') < 

ell, it 

< KaS^ \\{u,v,d) — {u ,v' ,d')\\^ 

^adj,±('W, V, d) - ^adj,±('w', d') < 

adj,± 

< Ka \\{u,v,d) - {u',v',d!)\\ . 



for m 



± 



Thus, since < (5 ^ cr, 



Q{w2) - Q{wi) 



< 2Ka\\w2 — wi I 



and therefore, for a small enough, it is contractive. 

The previous two propositions show that the operator Q is contractive. Let us denote by 
(u* ,v* ,d*) its unique fixed point in the ball B{2ki) C y. Now, it only remains to obtain the 
estimates stated in Lemma 6.4 The estimates for the hyperbolic variables are obtained as in 
the proof of Lemma |5.2[ For the elliptic ones it is enough to take into account that 



4 = Ck{Tj) = dk{Tj)e''^^ 

= geii,c,(0)(r,)e^^^ + {g,n^,^{u*,v*,d*){Tj) - geii,c,(0)(T,)) e^^^' 
= c^^ + {g,n,c,{u*,v*,d*){T,) - geii,c,(0)(T,))e^^^ 



and bound the second term using the Lipschitz constant obtained in Proposition 6.8 



We finish the section by proving Proposition 6.8, which completes the proof of Lemma 6.4 
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Proof of Proposition \6.8[ As we have done in the proof of Proposition 5.5 first, we stabhsh 
bounds for any {u,v,d) € B{2ki) C 3^ in the supremmum norm, which will be used to bound 
the Lipschitz constant of each component of Q. Indeed, if {u,v,d) G B{2ni) C 3^, it satisfies 



(89) and 



141 < KJ^^-'^/'' for k e Vf 

We bound the Lipschitz constant for each component of ^eii- We split each component of the 
operator between the elliptic, hyperbolic and mixed part. We deal first with the elliptic part. 
It can be seen that for k G V^, 

l^ciLc, {d'e'') - Z,n,c, {de'')\ <K^5^-'N{dk - 4) 

+ KJ Yl (4-4)- 



Therefore, 



Proceeding analogously, one can see also that 



< K^6^-''-NTj\\{u,v,d) - {u',v',d')\U. 



ell,± 



■j±2 



adj,± 



< K^S'-''NTj\\iu,v,d) - {u',v',d' 



Now we bound the mix ed t erms. Proceeding analogously and considering the properties of 



-Rmix.Cfc stated in Lemma 



6.2 



we can see that for m 7^ j it 2 
de 



ell,± 



1=1,2 



+ K^d^^^Sln'il/d) \\dk - 4||eii,± + KJ^^-'^'' Y - 4 



'^llell,± 



< K^C^^^ 5 ln\l/ 5) (1 + K,iV5(i-")/2^ \\{u,v,d) - {u,v',d')\ 
Therefore, using that 6 = e-^^ and (pi). 



Rnn.,c, {ue^\ve-^\ de'') - Rmi.,c, Ke^^*, v'e-^'-"', d'e'') ) dt 



ell,± 



< K„C'^^'^5\n^{l/5)\\{u,v,d) - {u,v',d') 
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So, we can conclude that for m G 



\\Geii,Ckiu,v,d) - gen,ck{u',v',d')\\^^^^^ < Ka5'^ ''In!^ {I /6)\\{u,v,d) - {u ,v' , d')\\*. 

Proceeding analogously we can bound the Lipschitz constant for Gen,cj±2 ■ bound it for m = 
j — 2, the other case can be done analogously. Here K denotes a generic constant independent 
of a. Note that now there is an additional term in i?mix.c _2- This implies that 



mix,Cj_2 ■ 



-J 

1/2 



-V3t 



i=l,2 



ladj,- 



1 /2 



< 



KaK'U^^^ (0(^)5)'^%-^* \\{u,v,d) - {u\v',d')\l . 



Therefore, integrating and applying norms, we obtain 



e"^* ( R 



mix,Cj_2 



ue-^\ve~^\de''] -R 



''mix,Cj„2 



d'e'' dt 



adj,- 

< Ka\\{u,v,d) - 



which leads to 



||^ell,c,_2(^i,'y,'^) - ^ell,c,_2(^t', v','^')||adj - - Ka\\{u,V,d) - {u ,v' ,d!)\\^. 

Now we bound the Lipschitz constant for the hyperbolic components of the operator. Note that 
we only need to bound the terms involving -Rmix.z since the other terms of the operator have 
been bounded in Proposition 5.5 We start with the Lipschitz constants of ^hyp,?;;- To this end 
we bound 



Rmi.,y^{ue^\ve-'^\d(^') - Rrai^,y^{ue^''\ve-^''',de''))dt 



< 



^* In,- - n' 



+ Yl + ^2))|4 - d'kldt, 



Ei'^'^i' Ei^^-^^i r* 



where we abuse notation concerning the 0{-) as before. Thus, integrating the exponentials and 
applying norms, one can easily see that 



Rmi.,yAue^\ye-^\de'') - R^^,y^{ue^''\ve-^'-'\de'')] dt 



V3t „,„-V3t ^jjt- 



< K^N 6'-' /6)\\{u,v,d) - {u',v',d')\U. 
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Therefore, applying norms and using condition on 5 from Theorem [Sj we obtain 



R 



mix,j/i 



ue^ ,ve 



de'' - R 



mix,j/i 



de" dt 



< K J'-' \n'{l /5)\\{u,v,d) - {u'v',d')\U 



3^^*, t;e-^*, de^*) - R^^,y^ (ne^*, ?;e-^*, de** 

Then, taking into account the results of Lemma |5.5[ one can conclude that 

< 



hyp,i;i 



hyp,-!;2 



hyp,t)i 



^hyp,Di (W, Vi d) - Ghyp,vi {U ,V ,d) 

< {(c^^H^^'^\n{l/5) + 5^-'\n^{l/5)^ \\{u,v,d) - {u ,v' ,d')\ 



Q\iyr>,V2 ^! d) - ^hyp,i,2 (""'> v\ d!) 



< 



hyp, 112 



< ^''^ ln(l lb) + (5i/2-2r in(l/5)^ d) - (u', f', d') |L 
Proceeding in the same way, one can obtain that 

^hyp,«i {U, V, d) - ^hyp,«i v, d') < 

hyp, til 

< (^{d^^'>6y^\n{l/6) + 6^-' ln\l/5)^ \\{u,v,d) - {u' ,v' ,d')\\^ 

Qhyp,U2 (u, V, d) - Ghyp,U2 {u , v' , d') 



hyp,U2 



< ln(l/5) + 5i/2-2r- ln\l/6)^ \\{u, v, d) - (n', v' , d')\\^ 



This completes the proof. 



□ 



7 The global map: proof of Lemma 4.8 



We devote this section to prove Lemma |4.8[ The continuous dependence with respect to initial 



conditions of ordinary differential equations gives for free that the map -S^j^^^ , defined in ( 66 ) , is 
well defined for points close enough to the heteroclinic connection defined in (33). Nevertheless, 



to prove Lemma |4.8[ we need more accurate estimates. 

7 

glob 

""loc 



Recall that the map ^S^jo^ defined in which is contained in M.(h) = 1 (see (31)). So, 



as we have done for B-} . we use the system of coordinates defined in Section 



4.1 



Recall that 



the initial section defined in (64), and the final section S^^^^, defined in ( |56[ ), are expressed 
in the variables adapted to the j*'* and (j + 1)** saddles respectively. Namely, in the coordinates 
{pf^ , , , qi'^ , c(^) ) and , ^^^^^ , 4^'+^^ , gf^^^ , c^^^^^ ) (see Section Q . To simphfy the 

exposition, first we will study the map iSgj^j^ expressing both the domain and the image in the 

variables {p^i \ (ii\p2 \ j c^-'-'). Then we will express the image of BL^^ in the new variables. 
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To simplify notation we denote the variables adapted to the j*'* and {j + 1)** saddles by 

{pi,qi,P2,q2,c) 

and 



Jj) Aj) U) (j) Jj) 

Pi tHI :P2 112 ' 



and we denote by the change of coordinates that relates them, namely 
Lemma 7.1. The change of coordinates is given by 



0^^(^1,91,^2,92,1 



UJP2 + UJ^q2 



-Ck 



for k e Pf+i U {j + 3} 



©i APi,qi,P2,q2,c) 

.7 — 1 



UJP2 + W 92 /, 2 



{uj Pi + uqi) 



0~,(Pi,9i,P2,92,c) 



-92 



0~ (Pl,9l,P2,92,c) = 



&~^{Pi,qi,P2,q2,c) = Rez+ — Imz 



V3^ 



&l^{Pi,qi,P2,q2,c) = Rez- — Imz, 



where uj = e^'^*/^ and 



r^ =1 



fcy^i-ljj + l 



(Pi + 9i - Pi9i) - {pI + 92 - P292) 



(98) 



^ =P2 + 92 - P292 



_Cj+2 



{u}p2 + lJ^q2) ■ 



Proof. We consider a point (p, g, c) and we express it in the new variables. We have to undo 
the changes (38) and (35) referred to the saddle j and then apply them again but referred to 
the saddle j + 1. The point (p, q, c) has associated variables r (as defined in (98)) and 6. We do 
not need to know the value of 6 to deduce the form of the change . Indeed, note that if we 
consider the changes (35) and (38) for the mode we have 



re 



id 



bj+i = Cj+ie = [uj P2 + ujq2 ) e 



id 



which implies 



UJP2 + 92 



(99) 



Using this formula and recalling that CfcC*^ = &fc = c^e*^, it is straightforward to deduce the form 
of Q~ for k G V^,^ U I? + 3). To deduce the form of Q~ and Q~ it is enough to consider the 
changes ( 35 ) and ( 38 ) for the mode bj to obtain 



re 



bj = cje 



ie 



[UJ 



''pi + ujqi) e'' 



Then, it is enough to use formula (99) to obtain and Gi . The others components can be 
obtained proceeding in the same way. □ 



55 



The next step of the proof of Lemma 4.8 is to express the section ^Y+i variables 

This is done in the next corohary, 



(pi, qi,P2,Q2, c) using the change obtained in Lemma 



7.1 



which is a straightforward consequence of Lemma 7.1 
Corollary 7.2. Fix a > and define the set 



where ^^+1 section defined in ( |56[ ) and 

yVj+i = < 1], \qi\ < r], \q2\ < rj, \ck\ < r] for k £ Vf and k = j ±2^ , 
Then, for r] > small enough, V^j+i can be expressed as a graph as 

P2 = w{pi,qi,q2,c). 
Moreover, there exist constants n' , k," independent oft] satisfying 

0<K < \/\-a'^ <K <\ 
such that, for any {pi,qi,q2,c) G VVj+i, the function w satisfies 

K < w{pi,qi,q2,c) < k". 
Once we have defined the section S^^^^, we can define the map 



*-'glob 



{pi,qi,q2,c) ^ ^giob(pi'9i'92,c) 



as 



*-'glob 



We want upper bounds independent of 6 and j for the transition time of the correponding orbits 
for this map. In the variables {pi,qi,P2,q2,c) the heteroclinic connection (33) is simply given 

by 

{p1it),q'i{t),P^2it),q!l{t),c\t)) = (^0,0,^^-^^=^,0,0) (100) 
(see jCKS+10] V Taking to such that 



1 + e2v^*o 



one can easily see that ^2(2^0) = \/l — o"^ and 2to ~ ln(l/a"). In the new coordinates this 
point is {pi,qi,P2,q2,c) = (0, a, 0, 0, 0) and thus belongs to the section qi = a. Then, thanks to 
Corollary 7.2, one can easily deduce that the time Tsj = Tsj (qi,Pi,P2, c) spent by the map 

''glob 



glob 



•^giob point iqi,Pi,P2,c) € Uj C is also independent of 5 and j. Recall that the 

difference between iSgj^j^ and B-^^^^^ is just a change of coordinates and therefore the time T^, 

spent by .Sgj^^jj is the same as Tgj . Thus, from know on we will only refer to T^j 





glob 



glob 



glob 



Next step is to study the behavior of the map B-'^^^^. In particular, we want to know the 
properties of the image set S'^^iohi^j)- 
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Proposition 7.3. Let us consider a parameter set Ij (as defined in Definition 4-6) and a I. 



product-like setUj. Then, there exists a constant K^j independent of j, N and 6 and a constant 
D^^^ satisfying 

such that the set B'l^^^{Uj) C satisfies the following conditions: 



CI 



SAT 



■^giob(^i) C D] X . . . X X Sj X D^+^ x...xB^j 



where 



{\ck\ < [U'il^ + Ky) <5(i-'-)/2} for k G 



C < \cj±2\ < 



1/2 



and 



1 /9 

{iPi,qi,P2,q2) G : \p,\, \q,\ < K^M^^l {c^^^s) , 
P2 = a, -D^^ 5 (ln(l/(5) - i?,) < 4^'^ < -D^^^ 5 (ln(l/<5) + }, 



C2 Let us define the projection ■k{p, q, c) = {p2, q2,Cj-2, ■ ■ ■ , cat). Then, 



V+2 



XO; X...XJ 



D(^)5(ln(l/5) - J(ln(l/J) + 1/K^) x{a}: 

where 

D^^__ = llci-')! < (m({; - ^^5^') <5(i-')/2} for k G 

1/2 



3- 



The proof of this proposition is postponed to Section 7.1 

Once we know the properties of the set &gio\^{l^j), there only remain two final steps. First 
to deduce analogous properties for the set S^gio^ji^^j) C S^^^. Second, to obtain a parameter set 
Ij+i and Xj+i -product like set Vj C S^^^^ which satisfies condition (70). These two last steps 



are summarized in the next lemma. Lemma 4.8 follows easily from it. 



Lemma 7.4. Let us consider a parameter set Xj+i whose constants satisfy 

/2 < C(^'+i) < 2D^^^ 
< mgl'^ < mjjp 
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and 



Then, the set 



M. 



(i+i) 



M. 



oU- 

(i+1) 



m 



M. 



ell,+ 
ell 

(i+1) 



adj, 



M. 



(i+1) 



adj,— 

(i+1) 

adj 

(i+1) 
hyp 



max <j M'i^_ + , kMX- } 



M. 



U) 



m 



m 



(j) 
ell 

(j) 
ell,- 



U) 
hyp 

max 



Vi+i = ^Lb(^i) n (P2, <?2, a, 5) = 0} , 



where gXj+i is the function defined in (94), is aXjj^i-product-like set and satisfies condition (70) 
Proof. It is enough to apply the change of coordinates given in Lemma 



7.1 



□ 



7.1 Proof of Proposition 7.3 



We spht the proof of Proposition 7.3 in several lemmas, which will give the needed estimates for 
the different modes. First, let us obtain rough bounds for all the variables, which will be used 



in the proofs of the forthcoming lemmas. Indeed, since we are restricted to M.{h) = 1 (see (31 )) 
we know that 

\Cm\ < 1. (101) 



Analogously, using the change ( 38 ) , one can see that 

\pi\ < 2, \qi\ < 2 for z = 1, 2. 
Now, we start by obtaining more accurate upper bounds for each mode. 



(102) 



Lemma 7.5. Consider the flow ^ associated to the vector field in (44) and a point {pi,qi,q2,cr,c) G 



Uj C S*-"*. Then, there exists a constant K„ > such that for t e [0,T„j ], <I>* (pi , gi , o", §2 , c) 

■' glob 

satisfies 

\^iM,qi,a,q2,c)\ < K„M^£j'^^-'-^/^ for m E Vf 

1/2 



<.±,(Pi,'Zi,^,92,c)| < kMJI± 



and 



|^pi(Pl,5l,f^,92,c)| < ^^<7^"hyp 



\^\, {Pi,qu^, q2, c) I < K^'I^X (^^'^'^ 



1/2 



1/2 



^L(Pl,9l,f^,92,c) -p^(t) 



€M^1u^,q2,c)\ < K^C^^^ 6 ln{l/ 5). 
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We defer the proof of this lemma to the end of the section. 



The bounds obtained in Lemma |7.5| are not enough to prove Proposition 7.3 since we need 
more accurate estimates for the ehiptic modes, the future adjacent modes and q2- We obtain 
them in the following three lemmas. 



Lemma 7.6. Consider the flow ^ associated to the vector field in (44) and a point {pi, qi, a, q2, c) G 



S'-"*. Then, there exists a constant Ku > such that for t G [0,T„j ] and k € , 



glob 



iT 



Proof. It is enough to point out that, using the bounds obtained in Lemma 7.5, the equation 



for Cfc in ( 44 ) can be written as 



where 7 satisfies ||7||oo < Kfj5^~'^~^^' . Then, to finish the proof of the lemma it is enough to 
apply the variation of constants formula and take into account that the time T„j has an upper 

''glob 

bound independent oi 5. □ 
Lemma 7.7. Fix values pi,qi,q2, cj -2 and Ck for k £ such that the set 

,i+2 



V = {ci,. . . ,Cj-2,Pi,qi,cr,q2} x D^. _ x {cj+3, . . .,Cj^} , 



where 



satisfies 



1,- 



\c,+2\<m^^l(c^'^S 



1/2 



VcUj. 



Consider the flow associated to the vector field in (44) and define the following map for points 
in T) 



T 



^adj ipi ,qi,cr,q2,c) = ^c.+l"'' (Pl ,Ql,cr,Q2, c) 
Then, there exists K^j > such that 

1/2 



Proof. Taking into account the estimates obtained in Lemma 7.5, the equation for Cj+2 in (44) 
can be written as 



_d / Cj+2 
dt \ Cj+2 



icj+2 - iuj {p^{t)Y Cj+2 + 7i+2(i) 
-iap^ + iu"^ {P2{t)) Cj+2 + lj+2{t) 



where P2 has been defined in (100) and 7 satisfies ||7llc>o < K^{C^^')d)^/'^5''' . Then, to finish the 
proof it is enough to apply the variation of constants formula. □ 



Now we obtain the refined estimates for q2- 
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Lemma 7.8. Fix values pi,qi,Cj±2 and cu for k G such that 



Q = {ci, . . . , c,_2,pi, gi, a}x [-C^^^6 (ln(l/5) - m^J^) , -C^^^6 (ln(l/5) 
satisfies 



+ m 



hyp 



x{cj+2,---,Cj^} 



Consider the flow associated to the vector field in (44) and define the following map for points 
in Q 



Fhypiq2) = iPi,Qi,cr, q2, c) 

Then, there exists K„ > and D^^^ satisfying 

such that 

-D^^h (ln(l/<5) - l/k^) , -D^^h (ln(l/<5) + l/K^)] C Fhyp(Q). 

Proof. Taking into account the estimates obtained in Lemma |7.5[ we write the equation for q2 
in as 

92 = Co(i)92 + Ci(i), 
where Co only depends on iii (100) and Qi satisfies 

llCllloo <i^.C(^')'^. 

Thus, the proof of the lemma follows from the variation of constants formula. □ 



We devote the rest of the section to prove Lemma 7.5 



J 

"glob 



Proof of Lemma 7.5. During the proof of this lemma the time t will always satisfy t G [0, 

and the norm || • Woo will always refer to the supremmum taken over this time interval . 

We start by obtaining the bounds for the non-neighbor elliptic modes. By (44), one can 
easily see that for k E 

d : |2_l/2 I 2 \ — 2 1 



dt 



Cfc_i + Cfc+i) Cfc - - (Cfc_i + Cfe+l ) Ck' 



Then, using (101), we have that 

'^1 |2 ^ I |2 

dt 

and therefore, applying Gronwall estimates we obtain that for t G [0,7^^^^^^], 

Proceeding analogously we deal with the adjacent elliptic mode Cj-2- Its associated equation is 



— I 

dr 



^2S--3Ci-2 + -Cj-3 Cj-2 

1/2 \ 2 2 ^ / 2 \ 2 2 

- X pi + uqi ) Cj-2 - - [ujpi +UJ qi) Cj_2 • 
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Taking into account the bounds in (101) and also (102), to obtain 

Cj-2\ < 5|Cj_2| 



5T , 

< e gi°b c 



■j-2\ 



dt 

which, applying Gronwall lemma, gives 

Analogously, one can obtain 

Now we obtain the bounds for the hyperbolic modes. We define 

Piif) = (^pi(Pi,9i,<^,92,c),$*^(pi,gi,cr,g2,c)). 

From (41), one can see that pi satisfies an equation of the form pi = Ai{t)pi where Ai{t) is a 
time dependent matrix (which of course depends on (pi , gi , o", ^2 , c) itself). Using (101) and 
(102), one can deduce that 

Pllloo < K^. 

Then, the fundamental matrix ^ satisfying ^'(0) = Id associated to this system satisfies ||^||oo < 
Kfj. Since pi can be just written 

pi(t) = ^(t)pi(O), 



using that by hypothesis |pi(0)|, |gi(0)| < M^^^ (C^^^d)"', we have that for t G [0,T„j ] 



hyp 



1/2 



glob 



1/2 



We finish the proof of the lemma obtaining the estimates for the {p2,Q2) components. To this 
end, let us point out that the equation for q2 can be written as 

q2 = ai{t)q2 + hi{t) 



where ai(t) and6i(t) are functions which depend on <I>p^ (pi , gi , cj, (72 , c) . Using (102) and the just 
obtained bounds for the non-neighbor and adjacent elliptic modes and for {pi-^qi) components, 
one can easily see that 

1 1 Oil 1 00 ^ and ||&i||oo 
Therefore, applying Gronwall lemma, we can deduce that 

\^\,{puqi.a,q2,c)\ < K ^C^^^ / 5) . 



To obtain the bounds for p2 we define ^ = P2 — P21 where P2 is the function defined in (100). 
Using (102) and (100) we have the a priori bound ||^||oo < 3. Therefore, from (44) we can deduce 
an equation for ^ of the form 

i = a2m + h2{t), 

where the functions 02 and 62 satisfy 

1 1 0-2 1 1 00 

< and II62II00 < K^5'^ . 
Then, applying Gronwall's lemma, we obtain 

llClloo < KJ"-' 

which implies the estimate for ^p^{pi, qi,a, (72, c) —^2- This finishes the proof of the lemma. □ 
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A Proof of Normal Form Theorem [2] 

To proof of Theorem [2j we consider as a change of variables F the time one map of a Hamiltonian 
vector field Xp, where F is the Hamiltonian 



ni,n2,n3,n4£Z2 



with 



Fnxn2n^rn — ] , ,n , i 1 19 rii — n2 + ^3 — 714 — 0, 

— \n2\ + |?T-3| — 1^4^ 

I i2 I i2ii i2 I i2/n 

\ni\ — \n2\ + \n^\ — \n/^\ 
-^nin2n3n4 =0 Otherwise. 

If we define the flow of the vector field associated to the Hamiltonian F, we have that 

UoT = H o^^^^_^ 

=n + {n, F}+ [ (i-t) {{n, f},f}o ^'pdt 

Jo 

=v + g + {v,F} 

+ {g,F}+ [\l-t){{n,F},F}o ^'pdt, 
Jo 

where {•, •} denotes the Poisson bracket with respect to the symplectic form = | J2n&z'^ ctn^a^- 
We define 



n = {g,F}+ [ {1- 1) {{n,F} , f} o 

Jo 



Then, it only remains to obtain the desired bounds for X-ji and T and to see that 

g + {V, F} = g. 

To obtain, this last equality, it is enough to use the definition for F to see that 

g + {'D,F}=^ ^ (l - «(|nip - |n2|^ + Insp - |n4|^)-Fni„2n3n4) aniO^anaO^ 

ni—n2+n3=n4 

=- y 

ni—n2+n3=n4 
|ni|2-|„2|2 + |„3|2 = |„4|2 



Now we obtain the bounds for Xn. We start by bounding X^g py, the vector field associated to 
the Hamiltonian {g,F}. We have to bound 
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Then, 



All the terms can be bounded analogously. As an example, we bound the first one, 



Y I^W^I <4 Y 

n,'in(il? n,m£Z 



E 



n\+n2=m-\-n 



E 



^^Y Y i«"iii««2i X] 



\(^n\ I |Ctn2 I I'-^ns I 



neZ2 ni+n2=n 
<0(ll«ll'i), 



meZ2 ni-n2+n3=m 



where, in the first line we have taken into account that \Fni 



n2n3n4 



< 1, and to obtain the last 
line we have used that each sum in the previous line is a convolution product. The other term 
in the reminder can be bounded analogously taking into account that 

^ (1 - t) {{n, F},F}o ^'pdt = ^ (1 - t) - g, f} o <^>'pdt 

+ [\i-t){{g,F},F}o^'pdt. 

Jo 

Analogously, one can obtain bounds for F — Id recalling that 

r = Id+ / Xpo^^pdt. 
Jo 

B Proof of Approximation Theorem [4] 

We devote this section to proof the Approximation Theorem |4| Throughout this section C 
denotes any positive constant independent of A'' and A. 

The solution /3'*' is expressed in rotating coordinates (see change (13)) and a is not. To 



compare them in a simpler way, we consider the equation (12) in rotating coordinates. To this 
end, we use that equation ( 10 ) also preserves the norm and therefore we perform the change 
of coordinates 

an = <7„e<«+l"l')*, (103) 
with G = — 2||a||^2- Then, the equation for g = {gn]n&^ reads 

-ign = £n{g) + Jn{g), 

where £ : ^ \s the function defined as 



(ni,?i2,"3)S,4(n) 



(104) 
(105) 



63 



with A{n) C (Z2)3 defined in and J : f £^ 

Hamiltonian 

where TZ is the Hamiltonian introduced in Theorem [2} Therefore, J satisfies 

\\Ji9)y=0{\\g\\l). 



is the vector field associated to the 



(106) 



Note that equation (104) and equation (14) only differ by J', that is, in the fifth degree terms of 
the equation. Moreover, note that g{0) = a(0) and therefore, by the hypotheses of Theorem |4| 

<7(0) = f3\0). (107) 

To prove that g and (3 are close we define the function ^ as 

^n=gn- Pn (108) 

and we apply refined Gronwall-like estimates to bound its £^ norm. Thanks to (107), we have 
that ^(0) = 0. Moreover, from equations (14) and (104), one can deduce the equation for It 
can be written as 



where 



Z^{t) =J 
Z\t) =DS (p^ 



Applying the (.^ norm to equation (109), we obtain 

d 



dt 



m\i^<\\z\t)\i, + \\z\t)i\\^, + \\z\i,t)\ 



(109) 

(110) 
(111) 
(112) 

(113) 



The next three lemmas give estimates for each term in the right hand side of this equation. 
Their proofs are deferred to the end of this appendix. 



Lemma B.l. The function Z^ defined in (110) satisfies < CX ^2 



Lemma B.2. The linear operator Z^{t) satisfies ||2^(t)^||^i < XlneZ^ fn{t)\^n\, where fn{t) are 
positive functions satisfying 

(114) 



fn{t)dt < CjN, 

where T is the time given in (22) and 7 is the constant given in Theorem^ 



To obtain estimates for Z'^{^,t) defined in (112), we apply bootstrap. 
Assume that for < t < T* we have 



£1 



< C\-^/^2~^. (115) 

A posteriori we will show that the time ( |22[ ) satisfies < T < T* and therefore the bootstrap 
assumption holds. 



64 



Lemma B.3. Assume that condition (115) is satisfied. Then the operator Z'^{^,t) satisfies 



Combining Lemmas B.l, B.2| |B.3[ equation (113) implies 
d 



dr 



To obtain bounds for we write this equation as 

d 
Jt 



5Af 



and we apply a Gronwall-like argument for each harmonic of Namely, we consider the following 
change of coordinates, 



(116) 



Then, we obtain 



J2 e/o(/"W+CA-V2),,^|^^| ^ ^^.5^ 



From this equation and taking into account that 

fn{t) + > 0, 

we obtain that 



5N 



dt 



Therefore, integrating this equation, taking into account that C(0) = ^(0) = and using the 
bound for T in (22) we obtain that 

To deduce from this bound, the corresponding bound for ||C||^i it is enough to use the change 
(116), the estimate (114) from B.2 and the definition of T in (22). Then, we obtain 

which implies 

Uy < 2e^^^||C||,i < 2e^^^A-325^7A^'- 

Therefore, using the condition on A from Theorem [4] with any k > C and taking N big enough, 
we obtain that for t G [0, T] 



and therefore we can drop the bootstrap assumption (115). 
Finally, taking into account ( |108[ ) and ( 103 ) we obtain 

E |«ne-<^+l"l')*-/3„| <CA-3/2, 

which is equivalent to statement ( 25 ) in Theorem |4j 
It only remains to prove Lemmas B.l B.2 and |B.3 
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Proof of Lemma B.l, Taking into account (106), we have that 

5 

£1 ■ 



rOI 



< c 



Therefore it only remains to obtain an upper bound for Taking into account that 

supp{/3^} C A, the definition of /J'^ in (23) and Theorem kil we have that 



N 



n.GA 



Now it only remains to point out that from the results obtained in Theorem 3-bis, we know 
that at each time all but three components of b are of size \bj\ < 5*^ for certain v > whereas 
the other two satisfy \bj\ < 1. Then, using the definition of 5 in Theorem [3| we obtain that 



N 



h (^"^0 \<Cil + Nd") < C, 



which implies 

This finishes the proof of the lemma. 



£1 



(117) 

□ 



Proof of Lemma B.2. To proof Lemma B.2 



this end, we use the function £ defined in ( 105 ) to obtain 



we start by analyzing each component of Z'^{t)^. To 



We define the functions /„, as 

fn{t) : 



(118) 



We analyze them differently whether n S A or n A. We start with the first case. 

We fix n G A and we want to study which terms in the right hand side of ( 118| ) are non 
zero. Indeed, each of the terms \d^^£k {f^^)\ is of the form Pnil^n2 with {iT'i,n'2,n) E A{k), 
{n,n2,ni) G A{k) ov ni = n2 = n = k (the last case arising due to the term — Ifl'nPfl'n in (105)). 
Then, these terms are non-zero provided ^ and f3^^ ^ 0. This condition is satisfied 
provided ni, n2 G A (see ([23])). Thus, we have that n, ni, n2 G A. Then, property 1a of the set A 
guarantees that A; G A. Properties 2a and 3a imply that n only belongs to two nuclear families. 
Therefore, it only interacts with seven vertices (recall that it can interact with itself through 
the term 



gn\ 9n in (105)). This implies that for a fixed n, 



except for seven values of fe, which correspond to the parents, children, spouse and sibling of 
n and n itself. Moreover, for the same reason, each term d^^E^ (/3^) which is non-zero, only 
contains a finite and independent of and n number of summands of the form /3ni/3n2 with 
(ni, 722, n) G A{k), (n, n2, ni) G A{k) or ni = n2 = n = k. 
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Reasoning in the same way, we can obtain analogous results for the terms |cfc £'„(/3^)|. 



From these facts, we can deduce formula (114) for n € A. Indeed, we have seen that /„ only 
involves seven harmonics of and that it is quadratic in them. Then, recalling the definition 
of in ([23l), Theorem 3-bis ensures that /n(t) has size /n ~ A ^ for a time interval of order 

- X-'^d" ~ X-^e-^"^ for the rest of 
X'^'yN'^. Therefore, 



A^ln(l/5) ~ X^jN (recall that -5 = e'^^) and has size /, 
the time, that is, for a time interval of order A^A^ln(l/(5) 



fn{t)dt <C{N + N^e-^"^) < C-fN. 

This finishes the proof for n E A. 

Now we need analgous results for n A. We need to see which terms of |9^„ffc {f^^) \ , that are 
of the form (3^-^(3^^, are non-zero. We know that they are non-zero provided (ni,n2,n) G A{k) 
or (n,n2,ni) G A{k) and ni,n2 S A. Note that know the case ni = n2 = n = k is excluded 
since n A and ni, n2 G A. Since n A and ni, n2 G A, property 1a implies that k ^ A. Then, 
property 6a guarantees that there are at most two rectangles with two vertices in A and two 
out of A. Therefore, we have that 

= 



except for three values of k, which correspond to n itself and the other vertex not belonging to 
A of each of these two rectangles. Reasoning as before each term d^^Sk [P^] which is non-zero, 
only contains a finite and independent of N and n number of summands of the form /Sn^/S^^ 
with ni,n2 G A. Then, reasoning as in the previous case, we obtain 



fn{t)dt < CjN. 



This finishes the proof of the lemma. 



□ 



Proof of Lemma B.3. To prove Lemma B.3, we split in \112\ as Z^ = Z{ -h Z^ with 



Zl{t) =£ (/3^ + i)-£ - D£ [p^) e 



Using the definition of £ in ( 105| ), it can be easily seen that 



,1 <C 



y ml + ml 



Then, using the bound for obtained in (117) and the bootstrap assumption (115), we 

obtain 



\zl\L<cx-'''^m\ 



£1 



We proceed analogously for Indeed, it satisfies 

5 



k=l 



and applying (117) and (115) again, we obtain 

II^IIL <CA-^/'lieil.i 
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Figure 4: Rectangles 



Thus, we can conclude that 



\z\,<c\-''i^m\ 



£1 



□ 



C A result for small initial Sobolev norm 

In Theo rem [T] we cannot ensure that the initial Sobolev norm ||u(0)||hs is arbitrarily small as is 
done in [CKS+lOj . One could impose this condition at the expense of obtaining a worse estimate 
for the time T. In this appendix we state an analog of Theorem [T] under assuming that ||ti(0)||//s 
is arbitrarily small. 

Theorem 7. Let s > 1. Then there exists c > with the following property: for any small 
jjL and large ^ ^ 1 there exists a a global solution u{t, x) of ^ and a time T satisfying 

cln(^//i) 

l^J 

such that 

\\u{T)\\h'^ > A and ||n(0)||//s < /i. 

Remark C.l. Combination of Theorems^ and^ covers all regimes studied in WKS+ldf . 

The proof of this theorem follows along the same lines as the proof of Theorem [T] explained 
in Section [3] taking /C = A/ fi. The only difference is the choice of the parameter A to ensure 



< T < 



\u 
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Indeed, as it is explained in Section [Sj we have that 



and, therefore, one needs to choose A such that A ~ /i. By Proposition 3.1 , the constant S; 



defined in (26), depends on A^. Nevertheless, in that theorem there is no quantitative estimate of 



this dependence. We will compute it here and show how it affects the estimates for the diffusion 
time T. 



We will show that there is a choice of the set A with S3 from ( 26 ) satisfying 



S3<B 



Ar2 



(119) 



for certain B > independent of N, e.g. B = QO^ applies. 

First, using this estimate we derive the time estimate in Theorem [T] from (119). Later we 
prove (119). We choose 



-B 



N2 



SO that A ^5*3 ~ fj,. Then, by Proposition 3.1 we have ~ In/C. Taking /C = A/fi, we know 



that there exists a constant c > such that 



A< 



and therefore, using formula (22) we obtain the estimate for the time. 



Now we prove (119). To this end we use the construction of the set A done in CKS^IO 



Recall that the authors first construct the set A inside the Gaussian rationals Q[i] and then 
multiplying by the least common multiple they map it to the Gaussian integers which is 
identified with Z^. Now, we want to place the points in Q[i] keeping track of the denominators. 
This gives us the size of the harmonics we are dealing with and, therefore, the size of 5*3. 

The placement of the modes in Q[i] is done inductively generation by generation. Namely, 
we first place Ai, then place A2 checking that the conditions 1a — 6a are satisfied, then place A3 
and so on. Note that the modes have to be close to the configuration called prototype embedding 
m ICKS+IO] . Sect. 4, since then we can ensure that (|15l) is satisfied. 

First generation: To place the first generation we consider a grid of points in <Q[i] with 
denominator 60^^. It is clear that we can place Ai in this grid with the points close to the first 
generation of the prototype embedding in [CKS"'"10] . It can be done so that (co)tangent of a 
slope between any two points in Ai has numerator and denominator bounded by Qi := 60^. 

Second generation: The set Ai is divided in pairs of modes which are the parents of different 
nuclear families. For each of these pairs, we need to place a pair of points of A2 forming a 
rectangle with the other pair. These new pair is going to be the children of the nuclear family. 
To place it we consider the circle C having as a diameter the segment between the considered 
pair in Ai. Then, the children have to be placed 



• at the endpoints of a different diameter of C. 

• they should belong to Q[i] and 

• the conditions 1a — 6a are satisfied. 
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Figure 5: Proper children's choice 



To see that the children belong to Q[i], we have to consider a diameter making a Pythagorean 
angle with the previous diameter, that is an angle 9 such that e*^ G (see Figure [5]). 

Let n = [\/ R/2\ be the integer part of \/ R/2. We lower bound the number of 0's whose 
tangent is rational with numerator and denominator bounded by R as ■\/R/2. To see that notice 
that any triple of the form a = m? — n^, h = 2mn, c = m? + with m < n is Pythagorean. 
Then there are n — 1 values for m giving a Pythagorean triple. 

Conditions 1a — 6a are satisfied provided the modes in A2 are not placed in certain points 
of the circle C. The number of these points is of order smaller than 60^. Indeed, we have to 
exclude: 

• The points of the previous generation (2^ points). 

• The points of A2 which have already been placed (at most 2^). 

• To avoid the existence of more rectangles besides the nuclear families, we proceed as 
follows. We consider 

— all the already placed points, 

— all the lines perpendicular to lines containing two of these points and passing through 
one of them, 

— all the circles having as a diameter the segment between two of the already placed 
points (see Figure [5]). 

Call C the set of these lines and C the set of these circles. The cardinality |£ U C| is at 
most of order 5'^. Then, we have to exclude all the intersections between any object in 
CuC with the circle C. 

• To ensure that condition 6a is satisfied, we consider the set V of the points which are the 
intersection between any two objects in >CUC It is easy to see that |P| is of order at most 
25^. Consider the sets 

— C containing the lines which are perpendicular to a line containing a point in V and 
an already placed point of A, and contain one of these two points. 
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— C containing the circles having as a diameter a segment whose endpoints are a point 
in V and an already placed point of A. 

The cardinality |£' U C'| is at most of order 60^. Then, we have to exclude also the 
intersections between an element in £ U C and C. 

We can place the children of the nuclear family at rational points of the circle C away from the 
ones just mentioned. To estimate its denominator we apply our estimate on the number of the 
Pythagorean triples. We have that the number of 0's with slopes whose tangent is given by a 
rational whose numerator and denominator is bounded by R lower bounded by \/R/2 — l. Thus, 
we can choose R = 60^^. Formula tan(Q + /3) = (tana + tan/3)/(l + tana tan /3) implies that 
(^2 < 2 60^^ Qi. Thus, denominators and numerators in Ai U A2 are upper bounded by Q2- This 
grid is accurate enough to place the pairs of A2 in the corresponding circles. Iteratively, we can 
place the following generations refining the grid at each step by dealing with Gaussian rationals 
whose (co)tangent has numerator and denominator bounded by 60^-^^ at the j generation. 
Therefore, after placing the generations and mapping the set A from Q[i] to Z[i] we obtain 
that all the modes n € A satisfy 

\n\<m^^\ 

This procedure can be done so that the final configuration of modes is close to the prototype 



embedding in ||CKS"'"10 to ensure that condition (15) is satisfied. Finally, to obtain the estimate 
(|119|), it is enough to take any B > 60^. 



D Notations 

/C — growth of the Sobolev norm of the solution from Theorem [T| 

s — index of the Sobolev space. 
T-L — the Hamiltonian of ([T]), defined in ([9]); 
T> — quadratic part of the Hamiltonian % defined in ([9]); 
Q — quartic part of the Hamiltonian % defined in ([9]); 

M. — abusing notation, mass of both the solutions of the equation ([T]) and of the toy 



model ( 18 ) 



{an(t)}„ez2 — Fourier coefficients of the solutions of ([T]) or, equivalently, solution of system 
On = 2ida:^n; 

r — normal form change for the Hamiltonian ([9]). It is given in Theorem [2] 
Q — resonant terms of Q. 

TZ — remainder (of degree 5) of the Hamiltonian T-L after performing one step of normal 
form, that is remainder of the Hamiltonian ^ o P. 

{an(t)}„gz2 — Solutions of the normalized Hamiltonian "H o P, given in Theorem [2| 



Ao{n) C (Z ) — collection of the resonance convolutions defined in (11); 
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{/3„(i)}ngz2 — rotated fourier coefficients, /3„ = a„e~*('-^"'"l"'l^)*. They satisfy (14). 
A{n) C (T?)"^ — collection of reduced resonance convolutions defined after (14); 

— 4 — number of energy cascades; 
A C essential Fourier coefficients given as a disjoint union of pairwise disjoint gen- 



erations: A = Ai U • • • Ajv- See Proposition 3.1 and preceding discussion 



solution to the Toy Model (18); 



h{b) — Hamiltonian of the Toy Model, given in (19); 



To — periodic orbits of the Toy Model (18) 



{c^^^lfc^j — coordinates adapted to the periodic orbit Tj after symplectic reduction, given 
m Section im 

{pi,qi,P2,Q2) — hyperbolic variables adapted to the periodic orbit Tj after diagonalization, 
given in Section [4. 1[ 



Zhyp^*, 2£^^, Zmix,* — types of remainder terms of the original Hamiltonian H after sym- 
plectic reduction and diagonalization near the periodic orbit Tj. Subscript means hyper- 



bolic, elliptic and mixed remainder respectively (see Lemma 4.1). 



S'jP — transversal section to the stable manifold of Tj, defined in (56) 



S° — transversal section to the unstable manifold of Tj, defined in (64) 



— map from S'P to S'^^i given by the flow of the Toy Model (18) (see Section 4 



Bf^^ — local map from to given by the flow of (18), defined in (65). 



"^giob — global map from S^^* to S^?^! given by the flow (18), defined in (66) 
a = 0{b) means |6| < Ka for some K independent of 6, a, N,j. 
a = Oa{h) means \b\ < Ka for some K independent of 5, N,j. 



^hyp — the change of coordinates for the hyperbolic toy model (see Lemma 5.1) 



^' — the change of coordinates for the full toy model (see Lemma 6.1). 

Rhyp,*, Rmix,*, , -2^^,* — Collection of remainder terms for the Full Toy Model after normal 
form transformation ^ (see Lemma 6.1). 

Vj C S^P — an open subset contained in the domain of deflnition of B^^^ so that B^^^iVj) C 



Uj C S?"* — an open subset contained in the domain of deflnition of 5f.,„v, so that 



■'glob'^''' 



glob 



Mj — initial conditions inside whose orbits under the flow have cancellation prop- 



erty (see Lemma 5.2) 
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yVj — an auxiliary set in the (p, g, c)-space (see Corollary 7.2) 



9l {P2,Q2,<7,S) — the cancellation function, defined in (94) and used in the definition of 



Nf 



To — time of evolution of the Toy Model in Theorem [Sj 
7 — constant which gives the relation between 5 and A^. 
K — constant from upper bound on time in Theorem [3j 



• A — rescaling parameter see (21); 



K — constant wich gives the relation between A and N . 



• T — time of evolution after rescaling, see ( 22 ) ; 



• {b^{t)}jLi rescaled solution to the Toy Model, given in (21); 



3A 



{/3„(t)}„g22 the lift of the above solution to the Toy Model to approximate solution to 



(14); 
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